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This  Is  the  final  technical  report  being  submitted  by  JAYCOR  to  the 
Defense  Advanced  Research  Projects  Agency  and  Office  of  Naval  Research 
under  Contract  N00014-84-C-0737*  "Concepts  for  Near  Continuous  Reception  of 
VHF  Signals  Using  Meteor  Burst  Propagation." 

This  program  has  been  sponsored  by  DARPA  through  the  Office  of  Naval 
Research  (Dr.  Fred  W.  Quelle*  Scientific  Officer).  The  total  cost  and 
fixed  fee  estimated  for  this  program  was  $240*057.  The  period  of 
performance  was  from  18  September  1984  to  17  May  1985.  The  final  report 
was  due  16  July  1985.  The  program  has  been  completed  within  budget  and  on 
schedule. 

The  objective  of  this  program*  as  given  In  the  original  JAYCOR 
proposal*  has  been  "to  design  experiments  leading  to  the  verification  of 
diversity  reception  techniques  for  tactical  C2  signal  collection."  The 
results  of  this  design  effort  are  summarized  In  this  final  report.  The 
experimental  program  that  has  been  designed  Is  referred  to  In  this  report 
as  the  VHF  Reception  Program  —  Phase  II.  Similarly,  the  program  that  has 
Just  been  completed  Is  referred  to  as  the  Phase  I  program. 

In  the  Phase  I  program,  we  have  Identified  three  key  physics  Issues 
concerning  diversity  reception  of  meteor  path  propagation.  These  Issues 
are:  ( 

•  //What  Is  the  gain  vs.  beamwldth  tradeoff  on  a  single  receiver 

antenna?_j 

•  1  Can  Independent  meteor  propagation  paths  be  obtained  at  a  single 

location  using  multiple  receiver  antennas?  < .  -* 

«l  What  spatial  separation  Is  needed  to  obtain  Independent  paths 
from  a  common  volume?  y  _ 

We  have  determined  In  our  Phase  I  program  that  these  questions  are  the 
principal  physics  Issues  remaining  to  be  resolved  before  a  netted  diversity 
reception  system  can  be  Intelligently  designed  for  testing  of  applications. 
The  Intended  applications*  which  may  Include  Interception  of  noncooperative 
VHF  communications  from  beyond  the  1 ine-of-slght*  will  not  be  directly 
addressed  by  this  program.  To  directly  address  applications  would  require 
a  program  designed  with  the  parameters  of  the  Intended  applications*  and 
would  require  a  much  greater  emphasis  on  signal  processing  and  data 
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analysis.  Such  a  program  would  attempt  to  demonstrate  diversity  gain 
(reliable  signal  from  Individually  unusable  signals)*  as  well  as  network 
gain  (simple  combining  of  signals  at  a  central  processor)  In  a  system 
closely  resembling  the  system  to  be  deployed  In  applications.  Instead*  the 
objective  of  Phase  II  experiments  Is  to  resolve  the  Important  physics 
Issues  with  the  simplest  experiments  and  supporting  modeling  that  can  give 
a  good  understanding  of  the  answers.  If  successful*  the  Phase  II  program 
would  provide  preparation  sufficient  to  begin  design  of  a  full-scale 
demonstration  of  an  efficient  VHF  diversity-reception  network  In  some 
subsequent  program. 

Consistent  with  our  perception  of  DARPA  objectives*  the  Phase  II 
program  Is  the  simplest  that  accomplishes  the  objectives  without  undue  or 
unacceptable  uncertainties  In  the  results.  The  Phase  II  program  features 
the  following  cost-saving  measures: 

•  A  single*  central  receiving  site  located  on  government  property 
(GFE  land  use)  within  ten  miles  of  the  site  where  most  of  the 
manned  test  operations  would  be  performed. 

•  Simple  unmanned  remote  transmitter  sites  located  on  government 
property  (GFE  land  use),  and  activated  by  telephone. 

•  Simple  signal  transmission,  signal  processing*  and  data  analysis, 
avoiding  a  costly  development  program. 

The  Phase  II  experiments  are  conceptually  simple.  To  address  the  gain 
vs.  beamwldth  tradeoff,  one  would  use  an  array  of  eight  receiving  antennas 
at  the  receiving  site  to  produce  three  distinct  receiving  beam  patterns 
having  widths  of  10°,  16°,  and  40°.  The  corresponding  antenna  gains  over 
Isotropic  are  18  dB,  16  dB,  and  12  dB,  respectively.  Percent-copy  would  be 
measured  for  each  beamwldth  pointed  at  several  locations  with  respect  to 
the  hot  spot  patterns.  To  address  the  question  of  Independent  meteor 
paths,  one  would  operate  several  distinguishable  transmitters 
simultaneously  at  distances  from  the  receiver  of  475  km  up  to  2000  km.  The 
transmitters  may  be  aimed  at  a  common  volume  or  at  separate  receiving 
sites.  The  correlation  of  the  transmitted  signals,  which  are  just 
distinguishable  tones  on  a  continuous  carriers  of  the  same  frequency  would 
provide  answers  In  the  simplest  manner  to  questions  of  spatial  and  angular 
diversity  and  path  Independence. 


The  Independence  of  meteor  paths  Is  an  unresolved  Issue  because 
experiments  to  date  have  typically  used  beamwldths  of  40°  to  50°.  With 
these  wide  beams  It  Is  known  that  even  at  distances  of  35  km,  antennas 
receive  virtually  Identical  Information.  The  Phase  II  program  would 
explore  the  possibility  that  narrow  beams  focussed  on  separate  hot  spots 
can  receive  uncorrelated  Information  from  independent  meteor  paths. 

The  gain  vs.  beamwldth  tradeoff  Is  unresolved  for  the  same  reason, 
namely,  that  narrow,  high-gain  receiving  beams  have  not  been  tried  In  VHF 
meteor  communications.  The  question  Is  whether  percent-copy  can  be 
enhanced  by  focusing  a  narrow  beam  onto  a  hot  spot.  Focusing  reduces  the 
solid  angle  subtended,  but  the  higher  gain  Increases  the  usable  number  of 
meteors  per  solid  angle.  In  Phase  II  one  would  experimentally  explore  this 
tradeoff  In  order  to  determine  the  optimal  receiving  antenna 
configurations. 

The  model  that  we  have  developed  In  Phase  I  (see  Section  2.2)  has  been 
exercised  to  simulate  the  Phase  II  experimental  conditions.  The  model 
presently  runs  on  an  IBM  PC  and  takes  5  to  8  hours  per  run.  The  runs  have 
been  done  overnight  at  no  direct  cost  to  the  Phase  I  program.  Because  the 
number  of  runs  performed  to  date  has  been  limited  by  time,  the  preliminary 
results  are  not  comprehensive.  The  results  of  the  code  runs  do  suggest, 
however,  that  one  may  face  some  counter-intuitive  surprises  In  Phase  II. 

At  first  we  expected  that  by  focusing  with  narrow  beams  on  separate 
hot  spots,  we  would  obtain  Independent  meteor  paths  from  each  hot  spot. 
The  code  suggests,  however,  that  there  Is  a  strong  correlation  between  the 
hot  spots,  apparently  owing  to  the  effects  of  the  antenna  sldelobes. 

We  expect  that  based  on  the  results  of  the  experiments  and  further 
modeling  In  Phase  II,  It  will  be  possible  for  communications  systems 
designers  to  design  the  optimal  receiving  antennas  and  system 
configurations  confidently  for  any  application  of  VHF  meteor 
communications. 

Section  2  presents  a  technical  discussion,  which  summarizes  much  of 
the  Phase  I  work.  The  model  developed  by  JAYCOR  Is  described.  Equipment 
requirements  are  delineated.  A  summary  of  our  research  on  data  correlation 
processing  Is  presented,  even  though  much  of  It  will  only  be  applied  In  a 
more  advanced  program  that  may  follow  a  Phase  II  program.  Section  3 
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discusses  the  technical  Issues  that  would  be  addressed  by  the  Phase  II 
program.  Section  4  describes  a  program  plan  for  Phase  II. 


2.0  TECHNICAL  DISCUSSION 


2.1  Meteor  Comaunlcatlons/Background 

As  meteors  enter  the  earth’s  atmosphere  they  form  long  columns  of 
Ionized  particles.  These  columns  diffuse  rapidly  and  usually  disappear 
after  a  few  seconds.  During  their  brief  existence#  however#  the  Ionized 
columns  will  reflect  radio  signals#  typically  for  hundreds  of  milliseconds# 
giving  rise  to  what  Is  variously  referred  to  as  meteor  scatter#  meteor 
propagation#  or  meteor  burst  communications.  The  radio  waves  can  be 
scattered  obliquely  to  provide  beyond  1  1  ne-of-s 1 g ht  VHF  radio 
communications  over  distances  of  as  much  as  2000  km  (VII  lard  et  al.»  1955). 
The  frequency  range  for  this  mode  of  communication  Is  30  to  100  MHz 
(Oettlng#  1980).  The  lower  limit  Is  the  frequency  at  which  the  effects  of 
Ionospheric  scatter  become  Important.  The  upper  limit  Is  set  by 
limitations  of  receiver  sensitivity. 

Meteor  burst  communications  has  a  number  of  applications  of  potential 
value  to  the  Department  of  Defense.  One  such  application  Is  the  reception 
of  noncooperative  tactical  VHF  communications  from  beyond  1 1  ne-of-s ight 
(BL0S)»  which  Is  Illustrated  schematically  In  Figure  2.1.  Depending  on  the 
percent-copy,  or  connectivity#  the  receiving  system  may  be  able  to  under¬ 
stand  the  message  or  context#  or  at  least  monitor  communications  traffic. 
Meteor  burst  communication  systems  have  also  been  considered  for  jam- 
resistant  and  Intercept-resistant  applications#  for  determining  the  loca¬ 
tion  of  noncooperative  transmitters#  and  for  emergency  communications 
networks.  In  this  last  regard#  an  Important  advantage  of  meteor  burst 
communications  Is  that  ’’the  nuclear  surv  1  v  1  ab  11  Ity  of  the  meteor  burst 
medium  Is  superior  to  other  BL0S  media  such  as  satellite  and  HF  skywave...” 
(Oettlng#  1980).  Ionospheric  disturbances  associated  with  nuclear 
explosions  often  greatly  Increase  D-layer  absorption,  which  reduces  the 
availability  of  HF  communications  relative  to  VHF. 

Research  In  meteor  burst  communications  has  been  sporadic  In  the  past. 
The  earl  lest  report  of  a  radio  effect  resulting  from  meteors  was  In  1931 
(Pickard#  1931).  However#  It  was  not  until  1946  that  a  clear  correlation 
between  radio  signals  and  Individual  meteor  trails  was  established  (Hey  and 
Stewart#  1946).  Numerous  field  experiments  on  meteor  burst  communications 
were  conducted  during  the  1950’s  (see  Bailey  et  al„  1955#  McKinley#  1961# 


METEOR  TRAIL  (80-120  km  ALTITUDE) 


Figure  2.1.  Schematic  of  YFH  Reception  Application 
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or  Sugar#  1964  for  brief  histories).  In  the  1960’s#  the  major  Investment 
In  long  range  communications  was  In  satellite  systems#  with  meteor  scatter 
development  being  largely  disregarded.  However#  the  advances  In  solid 
state  and  computer  technologies  In  the  1970's  resulted  In  meteor  burst 
communication  systems  becoming  competitive.  The  first  large  network,  the 
Department  of  Agriculture’s  SNOTEL  system#  was  developed  to  monitor  500 
remote  data  acquisition  stations  by  2  master  receiver  stations  using  meteor 
burst  communications  (Vancil#  1984).  Since  then#  numerous  commercial 
message  communication  systems  have  been  developed  using  meteor  scatter 
(NOSC,  1980). 

In  any  application#  an  Improved  percent-copy  confers  an  advantage  on 
the  receiver.  Near-continuous  reception  of  noncooperative  signals  is  not 
feasible  with  a  single  meteor-scatter  receiver#  because  meteor  burst 
communications  Involve  Inherently  weak-signal  systems.  Typically#  signal 
transmission  and  reception  Involves  90  dB  attenuation  of  signal  power#  with 
another  80  dB  or  so  lost  upon  scattering  off  a  meteor  trail. 

The  purpose  of  the  Phase  I  program  was  to  design  an  experiment  to  test 
systems  In  which  percent-copy  Is  Improved  by  diversity  reception. 
Diversity  reception  refers  to  any  technique  In  which  multiple 
antenna/receiver  paths  are  utl  1  Ized.  Simple  "network  gain"  is  achieved 
through  diversity  reception  by  combining  signals  from  all  the  receivers  at 
a  central  processor.  The  processor  gathers  the  received  signals 
Individually  and  "ors"  the  results.  In  this  simple  approach,  no  real 
diversity  gain  Is  achieved. 

In  Phase  I#  we  assessed  a  more  advanced  approach  to  signal  processing# 
In  which  "diversity  gain"  Is  achieved.  "Diversity  gain"  Is  the  extraction 
of  a  reliable  signal  from  Individually  unuseable  signals.  This  approach  is 
described  In  Section  2.4.  "Di verslty  gain"  wi  1  1  not  be  demonstrated  In 
Phase  II#  because  the  developmental  effort  required  In  data  correlation 
processing  was  not  deemed  to  be  warranted  until  the  significant  physics 
Issues  have  been  resolved. 

As  a  result  of  the  experimental  work  done  In  the  1950’s  and  early 
1960’s»  there  are  numerous  empirical  approximations  to  the  performance  of 
meteor  burst  communication  systems.  Two  approximations  for  the  performance 
of  systems  are  found  In  Oettlng  (1980)  and  NOSC  (1980).  Several  computer 
codes  are  In  use#  such  as  the  Department  of  Commerce  Meteor  Burst 
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Communications  Model  (Haaklnson#  1983)#  which  utilize  these  approximations 
to  estimate  the  throughput  of  a  system. 

The  empirical  approximations  to  the  performance  of  meteor  burst 
communication  systems  are  scaled  results  from  existing  experiments.  This 
approach  restricts  the  application  of  the  approximations  to  systems  similar 
to  the  experiments  used  as  the  basel  Ine  data#  and  from  which  the  results 
are  scaled.  Consequently#  only  simple  questions  relating  to  parameters 
affecting  the  performance  of  a  meteor  burst  communications  system  can  be 
answered  reliably  using  such  an  approach  to  modeling#  and  Inter¬ 
relationships  between  parameters  remain  uncertain. 

To  help  resolve  physics  Issues  and  design  effective  meteor  scatter 
experiments#  a  computer  code  was  developed  In  Phase  I  to  simulate  meteor 
burst  communication  links  realistically  and  reliably.  This  code#  and  the 
model  upon  which  it  Is  based#  Is  described  In  the  next  subsection.  The 
code  should  be  upgraded  and  extensively  Implemented  In  a  Phase  II  program. 
2.2  Modeling 

The  Meteor  Burst  Communications  Simulation  (MBCS)  code,  developed  by 
JAYCOR  In  Phase  I,  Incorporates  many  details  of  a  communications  system  to 
allow  the  Monte  Carlo  simulation  of  the  response  of  this  system  by  randomly 
generated  meteor  paths.  The  code  operates  on  an  IBM-PC  or  compatible 
microcomputer#  and  Is  self  prompting.  First#  the  user  Inputs  Into  a  file 
the  horizontal  beamwldth  and  offset  angle  for  both  the  receiving  and 
transmitting  antennas  and  the  vertical  beamwldth  and  takeoff  angle  for  both 
the  receiving  and  transmitting  antennas.  (The  code#  as  presently 
configured#  assumes  a  standard  dipole  antenna  radiation  pattern.)  Then# 
the  code  requests#  In  order: 

•  the  latitude  and  longitude  of  the  transmitter 

•  the  latitude  and  longitude  of  the  receiver 

•  the  Greenwich  Mean  Time  for  which  the  simulation  is  to  be 
performed 

•  the  date  for  which  the  simulation  Is  to  be  performed  (only  the 
month  and  day  of  month  are  required) 

•  the  frequency  (In  MHz)  at  which  the  transmission  is  to  be 
simulated 

•  the  threshold  at  which  the  receiving  amplifier  will  detect  a 
signal  (In  Watts) . 
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le  computer  code  opens  an  ascii  file  labeled  "mbrec.dat"  which  at  the 
jmpletlon  of  execution  contains  records  of  100  valid  meteor  reflections. 
ie  record  Includes  time  of  reflection*  peak  power*  and  duration  of  the 
Ignal  (time  the  signal  was  above  the  threshold).  The  code  also  simulates 
he  arrival  of  meteors  In  the  earth's  atmosphere  In  a  realistic  manner 
hlch  reflects  the  effects  of  the  day  of  the  year,  the  time  of  day* 
atitude,  etc.,  on  the  meteor  arrival  process.  The  output  then  Is  a  tlme- 
omaln  weighting  that  can  be  used,  for  any  given  input  signal,  to  compute 
he  received  signal.  With  such  a  computer  code,  detailed  questions 
elating  to  antenna  design  and  orientation,  diurnal  effects,  redundancy, 
andwldth,  etc.,  and  their  Inter-relationships,  can  be  studied. 

MBCS  was  developed  to  be  transportable  and  readable.  The  code  was 
eveloped  In  FORTRAN-77  on  an  IBM  PC  compatible  computer.  The  only 
xceptlon  to  the  use  of  FORTRAN-77  was  the  random  number  generator,  which 
ad  to  be  developed  In  Assembler.  However,  the  main  code  and  the  random 
umber  generator  will  operate  on  any  reliable  IBM  PC  compatible  computer, 
he  code  was  further  developed  following  the  precepts  of  structured 
rogrammlng  to  make  the  code  easily  understandable. 

This  subsection  describes  the  function  and  operation  of  MBCS  and  the 
eatures  of  the  underlying  model.  The  required  upgrading  of  the  code  and 
ts  Implementation  In  Phase  II  will  be  discussed  In  Section  3. 

Figure  2.2  Illustrates  the  geometry  used  In  the  model.  In  this 
:1gure*  the  following  lines  and  curves  are  illustrated:  Pj.  is  the  plane 
leflned  by  the  scattering  point  M,  the  transmitter  site  T,  and  the  receiver 
;1te  R.  The  geometrical  condition  of  scattering  is  that  the  meteor  trail, 
illustrated  In  the  figure  as  m,  be  In  a  plane  tangent  to  the  ellipsoid 
leflned  by  the  foci,  R  and  T,  and  the  point  M.  This  plane  Is  defined  as 
>2>  and  m  Is  a  ray  in  that  plane.  The  ellipse  Ei  is  defined  by  the 
Intersection  of  the  ellipsoid  and  the  plane  Pj_.  Clearly  Pi  is 
terpendlcu  1  ar  to  P2*  The  scattering  half  angle  is  measured  In  the  plane 
’j.  It  Is  Important  to  realize  that  the  geometric  condition  of  scatter 
Joes  not  demand  that  Pi  be  perpendicular  to  the  earth's  surface,  S.  If  P3 
Is  defined  as  the  plane  Including  the  verticals  at  R  and  T  through  R  and  T, 
:hen  P3  may  be  tilted  with  respect  to  P^  as  shown  In  Figure  2.2.  The 
>lanes  P^  and  P3  intersect  In  a  straight  line  through  R  and  T. 
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Figure  2.9.  Comparison  of  Simulation  Code  with  Experiments 
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jsed  rhombic  antennas  oriented  along  a  line  Intersecting  the  two  stations* 
with  their  main  lobes  Intersecting  at  a  point  100  km  above  the  earth  midway 
jetween  the  two  stations.  The  European  system  used  half-wave  dipoles*  with 
two  receiving  antennas  at  Winkfleld  so  that  phase  differences  could  be  used 
to  Infer  the  azimuthal  bearing  to  the  point  of  reflection.  Data  was 
:ollected  for  the  European  system  at  37  and  70  MHz. 

The  first  validation  run  studied  the  diurnal  cyclt.  effect.  Daily  data 
Is  published  for  the  NBS  system  In  Montgomery  and  Sugar  (1957)  for  the  time 
periods  6-7  May*  18-19  June*  and  11-12  August*  all  in  1952.  The  Eta 
Aquarld  meteor  shower  occurs  during  the  6-7  May  period,  and  the  Persld 
shower  occurs  during  the  11-12  August  period.  Since  no  unusual  meteor 
activity  occurs  during  the  18-19  June  period,  these  data  were  used  to  make 
comparisons.  The  NBS  data  were  collected  for  28  minutes*  and  extrapolated 
to  an  hourly  rate.  Assuming  a  Poisson  distribution  arrival  process  for  the 
meteors  during  the  28  minute  period.  It  Is  possible  to  compute  confidence 
limits  on  the  meteor  rate  published  In  Montgomery  and  Sugar  (1957).  Like¬ 
wise,  the  meteor  process  generated  by  the  computer  code  Is  simulated  as  a 
Poisson  process,  with  the  time  required  to  generate  100  valid  reflections 
(where  the  transmission  link  exceeded  a  given  threshold)  recorded,  and  so  a 
confidence  Interval  can  also  be  generated  for  the  computer  code  results. 
The  observed  and  generated  rates,  with  95%  confidence  Intervals,  are 
presented  In  Figure  2.9.  The  only  serious  difference  between  the  computer 
output  and  the  observed  data  Is  for  the  12:00  data.  As  the  geometries  for 
the  noon  and  midnight  periods  are  the  same,  comparable  rates  are  usually 
observed  for  these  two  periods,  which  tends  to  validate  the  computer 
results  and  to  suggest  a  possible  equipment  problem  for  the  observed  data. 

The  annual  cycle  In  meteor  rates  was  studied  next.  The  annual  cycle 
In  meteor  arrival  rates  most  typical  ly  quoted  is  that  of  Hawkins  (1956), 
which  separates  sporadic  rates  from  shower  rates.  However,  no  Information 
Is  given  about  how  the  data  were  collected,  so  It  Is  not  possible  to 
construct  confidence  Intervals  about  the  Hawkins  estimate  of  annual  cycle. 
Since  no  Information  Is  given  on  the  equipment  used  to  gather  the  Hawkins' 
data,  the  NBS  system  was  simulated  with  the  meteor  arrival  rate  adjusted  to 
give  the  same  annual  mean  number  of  meteors  as  reported  in  Figure  2.7. 
Figure  2.10  shows  the  meteor  rate  (dally  average  value  of  the  meteors  per 
hour)  for  mid-January,  March,  May,  July,  Sept.,  and  Nov.,  as  given  in 


12 


18 


00 


06 


12 


Figure  2.8. 


6  NOV  1956  7  NOV  1956 

(Fro*  Sugar,  Proc.  IEEE,  1964) 


Diurnal  Variation  of  Duration  of  Underdense  Meteor  Bursts 


22 


the  foci  R  and  T  (see  Figure  2.2).  A  valid  reflection  will  be  received 
when  the  angle  between  that  plane  and  the  meteor  path  Is  less  than 


0, 


m 


l(rt  4  rr> 

4RjRrCOS0 


.  2  2 

sin  0  cos 


0) 


where  L  Is  the  length  of  the  meteor  trail,  and  the  remaining  symbols  are  as 
defined  In  Figure  2.2  (Eshleman  and  Manning,  1954). 

The  Initial  radius  of  the  trail  Is  assumed  to  be  14  Ionic  mean  free 
paths  (Manning,  1958).  According  to  this  assumption,  the  initial  radius 
varies  exponentially  from  1.3  cm  at  80  km  to  12.6  m  at  120  km.  While  there 
Is  some  empirical  evidence  that  the  initial  radius  may  be  larger  at  the 
lower  altitudes  (Sugar,  1964),  the  distribution  of  signal  duration  was 
found  to  better  fit  the  data  using  the  Manning  hypothesis.  The  excellent 
agreement  of  the  MBCS  code  with  experiments  on  mean  duration,  which  Is 
about  to  be  demonstrated,  could  only  be  obtained  using  the  Manning 
hypothesis. 

Once  the  trail  Is  formed.  It  expands  by  diffusion  at  a  relatively  low 
rate,  producing  a  radial  distribution  of  material  that  Is  approximately 
Gaussian.  The  quantity  (4  Dt  +  rQ 2)1/2  may  5©  taken  as  the  approximate 
radius  of  the  trail  after  a  time  t.  Here  D  is  the  diffusion  coefficient 
and  rQ  is  the  Initial  radius  of  the  trail.  The  mean  dally  value  of  D 
varies  logarithmically  from  1  m2/sec  at  85  km  height  to  140  m2/sec  at  115 
km  (Greenhow  and  Neufeld,  1955).  As  the  variation  In  D  Is  In  fact  due  the 
variation  In  air  density,  there  Is  also  a  small  diurnal  and  annual  varia¬ 
tion  In  D,  which  has  not  yet  been  Incorporated  Into  the  computer  code.  The 
diurnal  variation  In  D  contributes  to  a  diurnal  variation  In  the  duration 
of  meteor  bursts,  which  Is  shown  In  Figure  2.8  (Sugar,  1964). 

The  JAYCOR  MBCS  code  has  been  benchmarked  by  performing  several 
simulations  to  compare  the  computer  code  results  with  published  data.  A 
difficulty  that  arose  was  finding  quality  data  combined  with  adequate 
system  details  to  perform  a  computer  simulation.  The  two  systems  that  were 
used,  both  of  which  involved  some  compromises,  were  the  National  Bureau  of 
Standards  system  operating  between  Cedar  Rapids,  IW,  and  Sterling,  VA 
(Montgomery  and  Sugar,  1957),  and  a  system  operating  between  Gibraltar  and 
Wlnkfleld,  England  (Bain,  1960).  The  NBS  system  operated  at  49.8  MHz.  It 
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coordinates  but  with  magnitude  42  km/sec»  and  rejecting  those  meteors  not 
moving  towards  the  earth. 

As  a  meteoric  particle  enters  the  earth’s  atmosphere  It  collides  with 
air  molecules#  which  then  become  trapped  In  Its  surface.  The  Impact  energy 
produces  heat#  which  evaporates  atoms  from  the  meteor.  The  atoms  move  off 
with  a  velocity  In  the  earth's  reference  frame  that  Is  substantially  equal 
to  that  of  the  meteor.  Collisions  between  these  high  velocity  atoms  and 
the  surrounding  air  result  In  the  production  of  heat#  light#  and 
Ionization#  distributed  In  the  form  of  a  long#  thin#  paraboloid  of 
revolution  with  the  meteoric  particle  at  the  head.  The  electron  line 
density  In  the  trail  Is  proportional  to  the  mass  of  the  particle  (see 
Figure  2.3  and  Table  2.1).  In  the  evaporation  process  each  Impact  of  an 
air  molecule  frees  many  meteoric  atoms.  Thus  the  total  mass  of  air 
molecules  striking  the  meteor  Is  small  compared  to  the  meteoric  mass.  As  a 
consequence  the  velocity  of  the  meteor  remains  quite  constant  unti  1  the 
meteor  Is  nearly  completely  evaporated  (Sugar#  1964). 

As  a  meteoric  particle  approaches  the  earth  no  appreciable  Ionization 
Is  formed  until  the  particle  enters  the  relatively  dense  air  at  heights 
below  about  120  km.  Above  that  height  collisions  of  the  particles  with  air 
molecules  are  not  frequent  enough  to  be  of  significance.  As  the  particle 
traverses  the  region  below  120  km#  It  vaporizes  rapidly  and  most  particles 
are  completely  evaporated  before  reaching  80  km.  The  higher  velocity 
particles  produce  trails  at  slightly  higher  altitudes.  Larger  mass 
particles  produce  trails  at  lower  altitudes.  Paths  with  greater  angle  of 
Incidence  also  tend  to  produce  trails  at  higher  altitudes  (Greenhow  and 
Hall,  1960). 

The  length  of  a  meteor  trail  Is  primarily  dependent  on  the  mass  of  the 
meteoric  particle  and  Its  angle  of  Incidence  onto  the  earth's  atmosphere 
(zenith  angle).  The  most  likely  length  for  a  trail  formed  by  a  sporadic 
meteor  Is  15  sec(z),  where  z  Is  the  zenith  angle  (Eshleman#  1957).  A  more 
realistic  distribution  of  lengths  will  be  Incorporated  Into  future  versions 
of  the  code. 

The  observability  of  the  forward  scattered  reflection  from  a  meteor 
trail  Is  a  function  of  whether  the  reflection  Is  within  the  principal 
Fresnel  zone.  The  center  of  the  principal  Fresnel  zone  occurs  where  a 
meteor  path  lies  exactly  within  a  plane  tangent  to  the  ellipsoid  defined  by 


The  computer  code  MBCS  uses  the  Inverse  power  distribution  to  randomly 
generate  particle  masses*  but  uses  the  relationship  between  mass  and 
density  given  In  Table  2.1. 

To  a  reasonable  approximation,  the  meteors  can  be  treated  as  being 
Isotropical  1  y  distributed  both  In  space  and  In  vel  oclty  to  slmul  ate  the 
meteor  arrival  process  at  a  fixed  time.  Although  there  Is  some  apparent 
variation  In  the  spatial  density  of  meteors  along  the  earth's  orbit,  as 
depicted  In  Figure  2.6,  this  Is  a  slowly  varying  annual  cycle  which  will 
not  change  over  any  smal  1  segment  of  time.  There  Is  also  evidence  of  a 
variation  in  the  density  of  the  meteors  about  the  ecliptic,  with  the 
meteoric  particles  tending  to  lie  In  the  same  ecliptic  plane  about  the  sun 
as  does  the  earth  (Bain,  1960,  Meeks  and  James,  1957).  The  former 
variation  Is  significant,  with  as  much  as  a  factor  of  2.3  Increase  In  the 
density  In  July.  Hence  this  factor  Is  accounted  for  In  the  computer  code 
as  affecting  the  rate  of  meteor  arrival  as  a  function  of  date.  The 
concentration  of  meteoric  particles  near  the  ecliptic  has  been  found  to  be 
a  minor  effect,  however,  with  the  Isotropic  model  giving  an  adequate  fit  to 
the  data  (Bain,  1960).  The  arrival  of  meteors  In  the  earth's  atmosphere  Is 
then  modeled  as  a  Poisson  arrival  process,  with  the  rate  of  arrival  a 
function  of  the  day  of  the  year.  The  number  of  valid  reflections  Is 
further  modified  by  other  pertinent  factors,  as  described  below. 

The  meteors  also  can  be  modeled  as  Isotropic  In  direction  of  travel; 
but  tend  to  have  velocities  near  42  km/sec,  which  is  the  escape  velocity 
for  a  particle  leaving  the  solar  system.  The  earth  travels  about  the  sun 
at  a  velocity  of  30  km/sec,  resulting  In  meteors  entering  the  atmosphere 
with  velocities  ranging  from  11.3  to  72  km/sec  (Sugar,  1964).  The  computer 
code  generates  random  velocity  vectors  uniformly  on  a  sphere  In  velocity 
space  of  radius  42  km/sec,  and  then  subtracts  the  earth's  velocity  vector 
to  obtain  a  geocentric  velocity  vector  for  the  meteoric  particle.  Those 
velocity  vectors  not  directed  towards  the  earth  are  then  rejected.  This 
method  of  considering  a  meteor's  velocity  relative  to  the  earth's  velocity, 
and  rejecting  all  meteors  not  moving  towards  the  earth,  properly  accounts 
for  the  diurnal  cycle,  latitude  effect,  and  seasonal  declination  effect. 
The  computer  code  accomplishes  this  by  generating  meteor  locations 
uniformly  on  a  sphere  In  conf Iguratlon  space  of  100  km  height  above  the 
earth's  surface,  generating  the  velocity  vectors  uniformly  In  polar 
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The  annual  cycle  Is  caused  by  two  factors:  the  earth’s  tilt  on  Its 
axis  and  a  variation  In  the  spatial  density  of  meteors  along  the  earth’s 
annual  path.  The  major  factor  Is  the  tilt  of  the  earth's  axis  relative  to 
the  ecliptic  plane*  and  Is  similar  to  the  comparable  factor  affecting  the 
diurnal  cycle.  The  relationship  of  a  location  to  the  earth's  ecliptic 
plane  Is  varied#  which  In  turn  affects  the  ratio  of  vertical  to 
horizontally  Inclined  trails  and  the  number  of  angles  of  Incidence  giving 
valid  reflections.  The  second  factor,  a  variation  in  the  spatial  density 
of  meteors#  Is  Illustrated  In  Figure  2.6  (Hawkins#  1956b).  Figure  2.7 
shows  the  annual  cycle  for  sporadic  and  shower  meteor  rates  (Hawkins# 

1956) . 

The  mass  distribution  of  sporadic  meteors  is  such  that  there  are 
approximately  equal  total  masses  of  each  size  of  particle.  The  approximate 
relation  between  particle  mass  and  dally  rate  Is  given  In  Table  2.1  (Sugar, 
1964#  VII  lard#  et  a!.#  1955)#  and  Is  111  ust  rated  In  Figure  2.3.  A  model 
for  the  probability  distribution  for  the  mass  of  meteoric  particles  which 
more  closely  agrees  with  the  data  Is  the  "Inverse  power"  distribution 
FM(m)  =  K  m-k  # 

where  k  =  2  for  sporadic  meteors*  and  k  =  1.5  for  shower  meteors  (Eshleman, 

1957) .  Note  that  Table  2.1  Is  equivalent  to  k  =  1. 

TABLE  2.1.  ORDER  OF  MAGNITUDE  ESTIMATED  METEOR  MASS  DISTRIBUTION 
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Figure  2.4.  Observed  Diurnal  Variation  of  Meteor  Burst  Rates 


Figure  2.5.  Predicted  Diurnal  Variations 


Meteors  can  be  divided  Into  two  classes,  shower  meteors  and  sporadic 
meteors.  Shower  meteors  are  col  lections  of  particles  al  1  moving  at  the 
same  velocity  In  fairly  well-defined  orbits  or  streams  around  the  sun. 
Their  orbits  Intersect  the  orbit  of  the  earth  at  a  specific  time  each  year. 
At  these  times  the  well  known  meteor  showers  are  observed.  The  shower 
meteors  account  for  only  a  small  fraction  of  all  meteors.  It  Is  the 
nonshower  or  sporadic  meteors  that  comprise  nearly  all  of  the  meteors  of 
Interest  In  radio  propagation.  Sporadic  meteors  do  not  have  well-defined 
streams,  but  seem  rather  to  move  In  random  orbits.  The  computer  code  only 
simulates  sporadic  meteors,  but  can  easily  be  modified  to  simulate  shower 
meteors  as  well. 

The  rate  of  Incidence  of  sporadic  meteors  at  the  earth  Is  further 
modified  by  two  cycles,  a  diurnal  and  an  annual  cycle.  The  first  of  these, 
the  diurnal  variation  In  meteor  arrival  rate.  Is  due  to  several  factors. 
The  major  factor  In  the  diurnal  cycle  Is  meteors  being  swept  up  by  the 
forward  motion  of  the  earth  during  the  morning  hours,  and  only  those 
meteors  overtaking  the  earth’s  motion  entering  the  atmosphere  during 
evening  hours.  This  results  In  a  maximum  occurrence  rate  around  6:00  a.m. 
and  a  minimum  rate  around  6:00  p.m.,  as  shown  by  typical  data  In  Figure 
2.4  (Montgomery  and  Sugar,  1957).  This  factor  Is  further  modified  by  the 
tilt  of  the  earth’s  axis  relative  to  the  ecliptic  plane,  as  this  varies  the 
location  that  Is  the  leading  point  on  the  sphere,  the  effect  of  which 
varies  with  latitude.  Ftnally,  the  orientation  of  the  receiver-transmitter 
path  determines  whether  the  horizontally  Inclined  meteor  paths,  which  tend 
to  occur  at  noon  and  midnight,  are  parallel  or  perpendicular  to  the  trans¬ 
mission  link.  Analytically  predicted  diurnal  variations,  accounting  for 
these  effects  are  presented  In  Figure  2.5  (Hines,  1956). 
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The  forward  scattering  cross  section  Is  that  of  a  cylinder  with 
effective  radius 
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The  transmission  equation  Is  then 
PR  GTGRX2sin2a 

PT  32  7r  2RtRr(  Rt+Rr)  ( 1-cos2  0  s  1  n2 0  ) 
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This  equation  applies  until  approximately  a  time 

2  ? 

rp  _  X  sec  0 

4  7T  20 

after  which  the  relation  for  underdense  trails  Is  applicable. 

It  Is  Interesting  to  note  at  this  point  a  violation  of  reciprocity  for 
transmission  In  both  underdense  and  overdense  cases.  The  power  received  In 
both  cases  Is  proportional  to  s1n2a»  where  a  Is  the  angle  between  the 
el ectr  1c  polarization  vector  at  the  meteor  trail  and  the  reflected  ray. 
This  term  appears  because  an  electron  reradiates  an  amplitude  pattern  of  a 
Hertzian  dipole*  which  McKinley  (McKinley,  1961)  describes  as  a  circle 
rotated  about  a  tangent,  or  a  holeless  doughnut.  The  axis,  along  which 
there  Is  no  scattered  field.  Is  coincident  with  the  electric  vector.  Upon 
Interchange  of  the  transmitter  and  receiver,  the  bilateral  symmetry  of  the 
transmission  equations  Is  In  general  broken  by  a  difference  In  the  angle 
a.  For  example,  for  vertical  ly  polarized  signals,  the  received  signal 
vanishes  for  meteor  bursts  directly  over  the  receiver,  but  has  maximum 
strength  for  bursts  over  the  transmitter.  The  orbital  motion  of  the  earth 
and  diurnal  variations  can  skew  the  number  of  meteor  Impacts  over  a 
transmitter  and  receiver,  particularly  on  east-west  paths,  so  that  the 
percent-copy  will  differ  If  receiver  and  transmitter  are  Interchanged. 
Preliminary  tests  using  MBCS  suggest  that,  upon  Interchange,  the  difference 
In  the  power  ratio  Pr/P.£  Is  typical  ly  of  order  10  to  20%,  but  can  be  much 
higher  under  certain  circumstances. 
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where 

X  is  the  wavelength# 

Ry  Is  the  distance  from  the  transmitter  to  M» 

rr  Is  the  distance  from  the  receiver  to  M, 

P  and  0  are  as  defined  In  Figure  2.2. 

As  the  trail  expands  by  diffusion#  the  phase  difference  of  the 
contributions  from  the  electrons  on  opposite  sides  of  the  trail  Increases# 
so  that  the  received  power  decreases. 

The  complete  transmission  equation  for  the  underdense  case  Is 
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where 


Py  and  PR(t)  are  respectively  the  transmitted  and  received  power# 

Gy  and  Gr  are  respectively  the  power  gains  of  the  transmitting  and 
receiving  antennas  relative  to  an  Isotropic  radiator  In  free  space# 
X  Is  the  wavelength  In  meters# 

re  =  2.8178  x  10-15m  Is  the  classical  radius  of  the  electron, 
r0  Is  the  Initial  trail  radius# 

q  Is  the  line  density  of  the  trail  In  electrons  per  meter# 
a  Is  the  angle  between  the  electric  polarization  vector  and  the 
reflected  ray# 

D  Is  the  diffusion  coefficient  In  square  meters  per  second#  and 
t  Is  time  measured  from  the  formation  of  the  trail  In  seconds. 

For  overdense  trails#  the  assumption  that  the  Incident  wave  passes 
through  the  trail  essentially  unmodified  Is  no  longer  valid.  For  this  case 
It  Is  assumed  that  the  wave  penetrates  the  trail  until  reaching  a  surface 
of  sufficiently  high  electron  density  to  be  reflected.  The  model  used  Is 
that  of  an  expanding  cylindrical  reflector  of  radius  rc.  After  a  time  the 
electron  density  everywhere  falls  below  the  critical  value  and  the 
underdense  model  Is  once  more  applicable.  However#  by  this  time  the  radius 
Is  quite  large  and  the  signal  contribution  is  rather  small. 
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The  distribution  of  energy  reflected  by  a  meteor  trail  Is  a  function 
Oi  many  variables.  The  Ionization  density  distribution  across  and  along 
the  trail*  the  orientation  of  the  trail*  the  VHF  wavelength*  the  polariza¬ 
tion  of  the  Incident  wave  relative  to  the  trail*  the  motion  of  the  trail 
either  as  part  of  the  process  of  formation  or  due  to  ionospheric  winds*  and 
the  straightness  of  the  trail  are  all  significant  (Sugar*  1964).  In 
discussing  the  reflection  properties  It  Is  convenient  to  divide  the  trails 
Into  two  classes*  underdense  trails  and  overdense  trails.  Underdense 
trails  are  those  In  which  the  electron  density  Is  low  enough  so  that  the 
Incident  wave  passes  through  the  trail*  and  for  which  the  trail  can  be 
considered  to  be  an  array  of  Independent  scatterers.  Overdense  trails  are 
those  In  which  the  electron  density  Is  high  enough  to  prevent  penetration 
of  the  Incident  wave.  Waves  reflect  from  overdense  trails  In  the  same 
sense  that  they  reflect  from  solid  conducting  cylinders.  At  the  long 
wavelengths  considered  by  this  code*  the  effective  duration  of  a  trail  Is 
long  compared  to  the  time  It  takes  the  trail  to  form,  and  the  trail  may  be 
considered  to  have  a  cylindrical  shape. 

The  mass  distribution  of  meteors  Incident  dally  on  the  earth’s 
atmosphere  Is  shown  In  Figure  2.3.  The  approximate  electron  line  density 
(In  electrons/meter)  produced  by  meteors  of  given  masses  Is  shown  by  the 
right-hand  scale  In  this  same  figure.  The  line  density  demarcating  under¬ 
dense  and  overdense  trails  Is  commonly  taken  to  be  q  10^  electrons  per 
meter  (Sugar,  2964),  as  shown  in  Figure  2.3.  The  code  treats  both  under¬ 
dense  and  overdense  trails. 

For  underdense  trails,  the  model  assumes  that  the  trail  is  an 
Infinitely  long  right-circular  cylinder  of  electrons  whose  diameter  Is  very 
small  compared  to  the  wavelength,  and  that  the  trail  electron  density  Is 
low  enough  that  the  incident  wave  passes  through  the  trail  without  major 
modification  (Sugar*  1964).  The  signal  received  can  be  computed  by  summing 
the  energy  backscattered  by  each  electron  In  the  trail  while  taking  proper 
account  of  the  phase  relations  of  these  contributions  (Lovell  and  Clegg, 
1948).  When  this  Is  done  It  Is  found  that  the  principal  contribution  of 
energy  from  a  trail  Is  from  Its  first  Fresnel  zone,  a  region  of  length 
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Hawkins  (Figure  2.7)  and  from  the  computer  code.  The  95%  confidence 

•  Interval  has  also  been  plotted  for  the  computer-generated  results.  There 
seems  to  be  acceptable  agreement  between  the  published  and  computed  annual 
cycles. 

The  peak  amplitude  distribution  for  several  reception  thresholds  (In 

•  microvolts)  Is  presented  for  the  NBS  system  by  Montgomery  and  Sugar  (1957). 
The  distribution  of  peak  amplitudes  generated  by  the  computer  code  Is 
compared  with  the  published  distribution  In  Figure  2.11.  Again,  generally 
good  agreement  Is  obtained  between  the  computer  code  and  published  data. 

•  The  duty  cycle  of  signal  transmission  (length  of  time  the  transmission 
link  remains  above  the  threshold)  Is  also  studied.  Montgomery  and  Sugar 
(1957)  present  the  distribution  of  time  the  transmission  link  Is  above  the 
threshold  for  the  NBS  system.  The  comparison  of  the  NBS  results  to  the 

4  computer  code  Is  presented  In  Figure  2.12.  Furthermore,  for  the  time 

period  used  In  the  study  (June  18-19),  the  NBS  system  had  a  mean  duration 
of  1.4  sec.  whl  1  e  the  computer  code  had  a  mean  duration  of  1.33  sec.  with  a 
standard  deviation  for  the  mean  of  .13  sec.  (This  agreement  provided  the 
q  rationale  for  using  the  Manning  hypothesis  for  radius  size,  as  opposed  to 

the  observed  data  reported  by  Sugar). 

A  critical  variable  for  some  applications  of  meteor  scatter 
communications  and  for  the  Phase  II  experiments  s  the  "hot  spot" 
£  phenomenon.  The  hot  spot  effect  Is  evinced  by  a  higher  density  of  valid 

meteor  reflections  off  of  the  line  between  transmitter  and  receiver  than  on 
the  line.  This  effect  Is  Illustrated  In  Figure  2.13.  Virtually  all  meteor 
paths  In  the  earth's  atmosphere  are  oriented  with  a  velocity  vector  pointed 
(  towards  earth.  A  reflection  occurs  at  that  point  on  a  meteor  trail  where 

the  path  to  the  meteor  trail  Is  orthogonal  to  the  lateral  component  to  the 
meteor  direction,  which  for  a  meteor  entering  the  earth’s  atmosphere  will 
only  be  on  the  center  line  If  the  meteor's  lateral  path  Is  exactly 
horizontal,  which  seldomly  occurs.  The  hot  spot  effect  makes  other  varl- 
|  ables  such  as  beamwldth  and  antenna  locations  (If  there  Is  more  than  one 

receiving  antenna),  as  well  as  diurnal  and  annual  variations,  very 
Important,  It  Is  also  Important  In  the  case  of  diversity  reception  systems 
(  with  several  Independent  receiving  antennas  oriented  towards  the  same 

transmitter  to  learn  the  effects  of  path  direction  and  diurnal  effects. 
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PERCENTAGE  OF  TIME  BURST  AMPLITUDE  EXCEEDS  ORDINATE 


PERCENTAGE  OF  TIME  BURST  AMPLITUDE  EXCEEDS  ORDINATE 


Figure  2.12.  Comparison  of  Simulation  Code  vlth  Experiments 

Duty  Cycle 


DISPLACEMENT  ALONG  GREAT  CIRCLE  (km) 

(Prom  Canad.J.Phys.  34) 

Figure  2.13.  Relative  Signal  Contributions  from  Various  Parts  of 

Meteor  Region  for  1000  km  Path 


[ 


Consequently*  the  detailed  predictions  of  the  hot  spot  effect  possible  with 
MBCS  Is  one  of  the  most  Important  capabilities  of  the  code. 

The  hot  spot  effect  Is  studied  In  detail  by  the  European  system 
described  In  Bain  (1960).  The  European  system  had  two  receiving  antennas 
separated  by  a  moderate  distance  to  allow  phase  differences  to  be  used  to 
estimate  the  bearing  to  the  point  of  meteor  reflection.  Very  definite 
diurnal  effects  were  found  In  the  hot  spot  locations.  Figure  2.14  shows 
the  published  and  computed  densities  of  bearings  to  meteor  reflections  for 
midnight*  6:00  a.m.*  noon,  and  6:00  p.m.  The  diurnal  hot  spot  effect  Is  a 
dynamic  one,  but  the  data  had  to  be  collected  over  an  Interval  of  time 
during  which  the  hot  spot  Is  changing.  Consequently*  the  distribution  of 
bearings  for  the  3  hour  periods  before  and  after  the  times  simulated  by  the 
computer  code  are  shown  In  Figure  2.14.  This  allows  one  to  develop  a 
picture  of  the  dynamic  effect  observed*  and  compare  this  with  the  effect 
computed  for  a  fixed  time  using  the  computer  code.  An  acceptable  agreement 
between  the  observed  and  computed  diurnal  hot  spot  effects  Is  seen  In 
Figure  2.14. 

The  above  comparisons  are  meant  to  be  a  benchmarking  rather  than  a 
definitive  validation  of  the  MBCS  code.  Much  Is  not  known  about  the  data 
sets  used.  A  more  comprehensive  comparison  Is  needed  using  data  collected 
with  known  systems.  However,  this  preliminary  comparison  between  published 
and  simulated  data  Is  certainly  encouraging,  and  suggests  that  the  Meteor 
Burst  Communications  Simulation  code  Is  a  promising  tool  for  the  design  of 
meteor  burst  communications  systems  and  experiments. 

2.3  EXPERIMENT  DESIGN 

2.3.1  Experiment  Objectives 

The  evaluation  of  the  VHF  diversity  reception  concept  has  brought  Into 
focus  basic  questions  concerning  VHF  meteor  path  reception  characteristics. 
The  objective  of  the  experimental  program  that  has  been  designed  is  to 
address  these  questions  relating  to  the  variation  of  reception  duty  cycle 
along  a  single  path  as  a  function  of  the  receiver  antenna  gain  and  beam- 
width*  and  the  relative  Independence  of  meteor  path  signals  as  a  function 
of  the  angular  and  spatial  separation  of  the  propagation  paths  in  multiple 
source  reception. 
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Figure  2.14.  Comparison  of  Hot  Spots  Simulation  Code  with  Experiments 


2.3.2  Experiment  Features 

The  experiment  will  provide  a  capability  to  evaluate  these  Issues  by 
allowing  a  variety  of  configurations  to  be  employed  using  realistically 
achievable  galn/beamwldth  receiving  antennas  at  a  centrally  located  site# 
and  (by  Invoking  reciprocity)  a  number  of  separated  and  remotely  controlled 
transmitters  at  a  variety  of  distances  and  angular  separations  from  the 
receiving  site.  (See  Figure  2.15). 

This  approach  allows  for  the  complexity  of  antenna  array 
reconfigurations  needed  to  vary  gain#  beamwldth  and  pointing  direction  of 
the  directive  antennas  to  be  done  at  the  central  manned  receiver  site.  The 
central  site  will  Include  the  needed  data  acquisition#  recording  and 
processing  analysis  efforts.  Transmitter  sites  will  be  relatively  simple 
and  can  employ  transmitters  controlled  via  a  commercial  dial-up  telephone 
line  to  each  location  (see  Figure  2.16  for  network  control  diagram). 

A  unique  feature  of  the  proposed  system  is  to  have  three  transmitters 
placed  at  some  combination  of  four  to  seven  licensed  locations  all 
operating  on  the  same  RF  carrier  frequency,  but  using  a  single-side-band# 
surpressed  carrier#  with  each  transmitter  radiating  a  distinct  audio  tone 
as  an  Identifier.  This  approach  allows  all  signals  to  be  received  on  one 
common  receiver  by  combining  signals  from  different  antennas.  Data  will  be 
recorded  as  audio  tones  (one  for  each  transmitters)  on  a  simple  analog 
recording  unit. 

2.3.3  Candidate  Site  Locations 

A  receiving  site  that  has  been  considered  Is  Camp  Elliot  near  Miramar 
Naval  Air  Station,  San  Diego.  This  area#  north  of  the  city  (within  ten 
miles  of  JAYCOR  headquarters  In  San  Diego)  Is  government  property  under 
the  control  of  the  Commander  NAS  Miramar.  A  request  for  use  of  space  and 
facilities  for  the  receiver  site  and  the  large  multi-element  antenna  arrays 
has  been  made.  A  response  from  the  Publ  1c  Works  Office  at  Miramar  gives 
preliminary  approval  of  the  request  pending  detailed  selection  of  a 
specific  site  and  clarification  of  environmental  and  operational 
constraints.  (See  attached  letters  at  end  of  Subsection  2.3#  JAYCOR  to  NAS 
Miramar  and  NAS  Miramar  to  JAYCOR.) 

Four  prospective  transmitter  site  locations  have  been  Identified  which 
have  the  features  of  being  appropriate  distances  from  the  receive  site  (400 
to  1400  km)  and  lying  on  different  angular  paths  (north-south  to  east-west) 


f0  T1  TN  T2 
TRANSMITTERS  ALL  AT  fQ  +  TONE 


LIVERM 


COLORADO 
\  ...»  SPRINGS 
A'  T, 


V  *£  /sf 

V.  .*o 

V /  ..*•• 

SAN  DIEGO®*4** - - 


V.****  v#  ^  V  albUerque 


— 1000  km 


ANTENNA  ARRAY  #1 


ANTENNA  ARRAY  #2 


RECEIVER  SITE  WITH  ANTENNA  FARM 
•3 ARRAYS  DIRECTED  AT 

1)  Tx 

2)  T- 


3)  T 


RE-08288  Figure  2.15.  VFF  Reception  Prograe  —  Phase  II» 

Exper leant  Configuration 


1200  Hz 


Direct-dial  telephone  to 
site  from  the  receiver  s 
deactivates  transmitters 


from  San  Diego.  These  sites  are  also  easily  accessible  to  facilitate 
repair  or  maintenance  service.  Additional  transmitter  sites  should  be 
licensed  to  be  available  for  the  temporary  placement  of  a  transportable 
transmitter  to  Increase  diversity  of  distance/angle  data.  The  four  primary 
candidates  transmitter  sites  (see  map  In  Figure  2.15)  and  their  ranges  from 
San  Diego  are: 


Camp  Parks  (Livermore#  CA)  800  km 

Nellis  AFB  (Las  Vegas,  NV)  475  km 

Fort  Carson  (Colorado  Springs#  CO)  1300  km* 

Klrtland  AFB  (Albuquerque#  NM)  1000  km 

Criteria  for  additional  transmitter  site  locations  are: 

1)  A  site  col  Inear  with  an  existing  site  but  at  different  range. 

2)  A  site  within  a  few  angular  degrees  of  an  existing  site. 

3)  A  site  to  evaluate  maximum  viable  range  (l.e.#  1500-2000  km). 

Prospective  transmitter  sites#  which  can  be  selected  from  existing 

Army/Navy/Air  Force  and  government  property#  will  require  a  minimum  of 
space  and  support  facilities.  An  example  of  the  present  facilities 
operated  by  the  Air  Force  Is  shown  In  Figure  2.17. 

2.3.4  Frequency  Selection  and  Allocation 

The  potential  application  of  the  proposed  concept  Is  for  the  lower  VHF 
band  from  about  30  to  100  MHz.  Below  30  MHz  there  Is  a  high  probability  of 
long-range  Ionospheric  propagation.  Fortunately#  the  solar  cycle  Is  at  a 
low  for  the  next  few  years#  and  should  minimize  low  VHF  Ionospheric  path 
signals.  The  desired  frequency  for  this  experiment  Is  In  the  30-50  MHz 
range.  Above  50  MHz,  the  duty  cycle  Is  much  lower  than  In  the  low  VHF 
range. 

This  low  VHF  band  Is  assigned  In  one  to  two  megahertz  segments  to  the 
U.S.  Government  and  non-government  agencies  and  Is  primarily  used  for 


*  A  letter  of  tentative  approval  from  Fort  Carson  is  attached  at  the  end  of 
Subsection  2.3. 


local -area  master-station  communications  to  multiple*  remote  (often  mobile) 
sites.  The  Federal  Government  frequency  assignments  are  administered  by 
the  Intradepartmental  Radio  Advisory  Committee  (IRAC)  (ref.  Telcon  with  Mr. 
John  Caufey  202/724-3346)  In  coordination  with  the  FCC.  Selected  segments 
of  the  low  VHF  band*  1.e.»  33-34  MHz*  37-38  MHz*  40-42  MHz  are  devoted  to 
military  tactical  radio  and  some  known  government  meteor  communication 
users  such  as  the  Department  of  Agriculture  SNOTEL  network*  the  Federal 
Emergency  Management  Agency  (FEMA)  civil  defense  system,  an  Alaska  Meteor 
Burst  Communications  System  for  federal  bureaus  during  auroral 
disturbances,  a  commercial  Northern  Natural  Gas  of  Omaha  Oil  Pipeline 
Communication  System  In  the  midwest*  and  a  Canadian  Meteor  Burst  Store  and 
Forward  System  (JANET)  Initiated  In  the  mid  50’s. 

A  frequency  search  of  the  non-government  Master  Frequency  List  was 
made  on  the  microfiche  list  at  the  San  Diego  field  office  of  the  FCC.  All 
frequencies  In  the  lower  VHF  band  appear  to  be  assigned  to  low-power 
unscheduled  users  such  as  city,  county,  and  state  governments  and  private 
users  for  network  communications  between  a  central  office  and  remote 
stations  and  a  variety  of  police*  sheriff*  civil  defense  and  fire 
department  communications  on  the  FM  channels  assigned.  A  single  frequency 
often  has  a  number  of  different  users.  The  government  Master  Frequency 
List  Is  maintained  at  area  regional  centers  and  Is  a  classified  document. 
2.3.5  Procedure  for  Authorization  of  U.S.  Governmental  Frequency 
Assignment 

Assignment  of  a  frequency  to  perform  a  government  sponsored  experiment 
Is  Initiated  by  submitting  a  government  form  DD  1494  from  the  sponsoring 
agency  through  an  area  FCC  office.  The  Information  required  on  this  form 
Includes  specifics  of  the  transmitter  locations,  transmitter  power  (ERP), 
antenna  types*  modulation  bandwidth*  operating  schedules*  etc.*  and  a  list 
of  the  equipment  types  to  be  used  In  generating  and  radiating  the  RF 
signals.  Discussions  with  the  Western  Area  Frequency  Coordinator,  Mr. 
Turklewcz  at  Pt.  Mugu  California  (805/989-7  983)  and  with  the  San  Diego  FCC 
Coordinator,  Mr.  Wayne  Smith,  Navy  Comm.  Station*  SD  (619/696-5621) 
Indicate  frequency  coordination  and  assignment  typically  take  from  90  to 
180  days,  and  frequency  assignment  Is  based  on  a  priority  system  assigned 
to  the  government  sponsoring  agency.  (See  sample  DD  1494  In  attachments  at 
end  of  2.3.)  Because  this  experiment  would  employ  multiple  transmitters 
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jperatlng  on  the  same  frequency  at  widely  dispersed  locations  (typically 
>ne  In  each  of  3  to  5  western  states)*  frequency  coordination  may  be  fairly 
;omplex,  and  assignment  processing  time  may  be  greater  than  normal. 
Authorizations  will  be  sought  for  the  prospective  locations  previously 
nentloned  and  for  alternate  locations  in  the  event  that  frequency  conflicts 
lrlse. 

2.3.6  Experiment  Plan  and  Configuration 

The  experiment  will  involve  collecting  data  periodically  to  measure 
joth  diurnal  and  seasonal  effects  of  three  different  VHF  meteor  path 
:on  f  1g  u  rat  1  on  s.  These  three  experiment  configurations  and  the 
:orresponding  test  objectives  and  Implementations  are  briefly  described 
below  and  refer  to  figures  and  sketches  which  describe  the  antenna 
features*  propagation  paths,  and  network  configurations. 

1.  Gain  vs.  Beamwldth  (Figure  2.18) 

Objective:  Determine  the  signal  reception  duty  cycle  as  a 
function  the  ga 1 n/beamw 1 dt h  of  a  receiving 
antenna  when  directed  towards  an  over-the-horizon 
signal  source. 

Approach:  Ga  1  n/Beamwl  dth  variations  (see  sketches  of 

antenna  configurations.  Figure  2.19). 

Forward  Gain  Horizontal  (azimuthal)  Beamwldth 
12  dB  40° 

16  dB  16° 

18  dB  10° 

Data  Base:  Record  detected  signal  received  during  specified 
period  of  time  (duty  cycle).  Data  taken  at  same 
times  each  day  from  two  different  transmitter 
stations  to  acquire  two  Independent  sets  of  data. 

2.  Angular  Diversity  (Figure  2.20) 

Objective:  Determine  the  dependence/Independence  of  signal 
reception  of  two  signals  in  the  same  geographic 
area  when  antennas  look  at  different  meteor  path 
volumes . 
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AEALY  TO 
ATTENTION  Or 


DEPARTMENT  OF  THE  ARMY 
headquarters,  FORT  CARSON 
AND 

HEADQUARTERS,  4TH  INFANTRY  DIVISION  (MECHANIZED) 
FORT  CARSON,  COLORADO  80913-5000 


June  6,  1985 


Facility  Management 


Mr.  George  F.  Ross 
Senior  Engineer 
JAYCOR 

Hilton  Parkway  *2 

A570  Hilton  Parkway 

Colorado  Springs,  Colorado  80907 

Dear  Mr.  Ross: 

With  reference  to  a  phonecall  on  June  A,  1985  between  yourself 
and  William  Davis,  Chief,  DEH  Facility  Management,  Fort  Carson  would 
entertain  a  proposal  to  temporarily  site  your  project  on  this  installa¬ 
tion.  In  all  probability  the  land  could  be  made  available,  but  an 
available  enclosed,  lockable  building  that  has  electrical  power  to 
it  may  prove  more  difficult  to  provide. 

I  would  ask  that  you  contact  Mr.  Davis  at  579-3038  to  set  up 
a  meeting  here  at  Fort  Carson  to  further  discuss  the  details  of: 
your  proposal  and  to  visit  potential  sites. 

Sincerely, 


HlkrWi  Bfcown,  Jt 
Colonejt,  Engineer 
Director  of  Engineering  and 
Housing 
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IN  REPLY  REFER  TO 


DEPARTMENT  OF  THE  NAVY 

NAVAL  AIR  STATION  MIRAMAR 
SAN  DIEGO  CA  92  145 

11011 

Ser  1832C/g  ^  j  g 

JUN  1  4  1985 


Mr.  Ralph  M.  Wheeler 
JAYCOR 

11011  Torreyana  Road 

P.0.  Box  85154 

San  Diego,  CA  92138-9259 

Dear  Mr.  Wheeler: 

We  are  in  receipt  of  your  letter  of  16  May  1985  requesting  space  and  limited 
facility  support  for  a  temporary  antenna  farm  located  on  East  Miramar  land. 
Your  request  has  been  reviewed  and  is  satisfactory  in  principle  to  this 
command,  however,  additonal  information  must  be  considered  before  final 
approval  can  be  granted. 

At  your  convenience,  1  would  like  to  request  a  meeting  between  your  staff  and 
mine  to  discuss  items  such  as  site  location,  utility  requirements, 
environmental  impacts,  aircraft  navigational  constraints,  security  and  license 
agreements.  Please  contact  either  Ms.  Merrily  Severance  or  Mr.  Rene  Trevino 
of  my  staff  at  (619)  271-3321,  to  arrange  for  this  meeting. 

As  standard  practice,  we  also  request  that  you  have  the  Office  of  Naval 
Research  contact  us  to  verify  your  government  contract. 

Sincerely, 


NAS-MiR. 5216/16  (6  621 
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AC  power:  120  volts  AC,  1000  watts  maximum. 

Operational  schedule:  2-3  hours  per  day,  5  days  per  week,  for  up  to  12  months. 

|  Security  classification:  UNCLASSIFIED. 

Physical  security:  Locked  room  for  receiver  equipment. 

Maximum  number  of  personnel  involved:  Normal  operation,  2  people;  installation  and 
set  up,  3  to  5  people. 

§  Program  proposed  under  government  contract  N00U14-84-C-0737. 

Office  of  Naval  Research 
Arlington,  VA  22217 
Contract  Monitor:  Mr.  Fred  Quelle 
•  Telephone:  (617)  451-3171 

Upon  your  approval,  we  would  select  the  preferred  building  and  site  and  would 
coordinate  our  test  activity  with  your  representative. 

g  A  response  to  this  request  and  any  questions  regarding  this  effort  should  be 

directed  to  me. 

Yours  truly, 

Ralph  M.  Wheeler 
VHF  Experiment  Director 
9  (619)  453-6580  ext.  220 

RM  W:sb 
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d/UJCOR 


May  16,  1985 


Capt.  Cary  E.  Hackanson 
Commanding  Officer 
Naval  Air  Station,  Miramar 
San  Diego,  CA  92145 

Dear  Sir: 

A  concept  to  evaluate  the  reception  of  VHF  (40  mHz)  signals  from  over-the-horizon 
is  being  investigated  by  JAYCOR  under  a  government  contract.  An  experiment  is 
planned  in  which  we  would  listen  (receive  only)  to  a  far-off  transmitter  (100's  of 
kilometers).  Special  antennas  would  be  used  which  require  a  space  about  the  size  of  a 
football  field.  Available  government  space  and  some  limited  facilities  support  (AC 
power  and  housing  for  radio  equipment)  was  suggested  by  our  government  contract 
monitor  in  order  to  keep  costs  of  the  program  to  a  minimum. 

At  the  suggestion  of  your  public  works  people,  an  examination  of  useable  areas  at 
Camp  Elliot  was  made.  The  following  facilities  and  space  appear  to  satisfy  our 
requirements: 

Building  E-137  and  the  open  land  to  the  north 
Building  E-133  and  the  open  land  to  the  north 
Building  E-131  and  the  open  land  to  the  south. 

We  request  your  approval  to  propose  the  use  of  one  of  these  locations  for  conduct  of  this 
experiment.  We  fee  I  confident  that  this  experiment  could  be  successfully  completed  on  a 
*no  interference*  to  NAS  activities.  The  following  requirements  are  stated  for  your 
evaluation  of  this  request: 

Receive  only  operation  (no  transmitted  RF  signals  would  be  radiated. 

Equipment  quantity:  1  or  2  receivers 
Receive  antenna  features: 

Array  of  2  to  6  antennas  with  maximum  length  50  feet  and  width  of  150  feet. 
Maximum  height  35  feet  (This  is  lower  than  the  local  power  lines  in  the 
requested  area). 

Inside  space  required:  200  to  300  square  feet  of  warehouse  or  office  space. 


11011  Torreyana  Road  —  Post  Office  Box  85154  •  San  Diego.  California  92138-9259  •  (619)453-6 


ATTACHMENTS  TO  SUBSECTION  2.3 


o  Letter  JAYCOR  to  NAS  Mframar 
o  Letter  NAS  Miramar  to  JAYCOR 
o  Letter  Fort  Carson  to  JAYCOR 
o  DD  Form  1494 


Figure  2.21.  Spatial  Diversity  Experiment 
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Approach:  Use  directive  antennas  to  point  at  different  hot 

spots.  Antenna  beamwldth  will  be  1  ess  than  1/2 
the  angular  separation  of  the  optimum  hot  spots. 
Two  receivers  (one  rented  for  2  months)  will  be 
used  In  this  experiment  with  an  offset  In 
frequency  to  differentiate  signals. 

Data  Base:  Time  related  comparison  of  received  signals  from 
each  source. 

3.  Spatial  Diversity  (Figure  2.21) 

Objective:  Determine  the  dependence/Independence  of  signal 
reception  of  signals  located  at  various  distances 
and  angular  separations  from  the  central  receiver 
location. 

Approach:  Use  narrow  beam  antennas  to  point  at  optimum  hot 

spots  on  each  link  from  sources  located  at 
different  angles  and  distances  from  the  reception 
point. 
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CONFIGURATION  1  ANTENNA  CONFIGURATION  2  ANTENNA  CONFIGURATION  3 


Figure  2.19.  Galn/Beamrldth  Antenna  Array  Configurations 


Figure  2.18.  Ga1n/Beamw1<fth  Experiment 
Two  data  sets  from  two  separate  sources 


APPLICATION  FOR  FREQUENCY 
ALLOCATION 


A.  TARGET  STARTING  OAT*  TOR  EACH  STAGE. 


CLASSIFICATION 

PAGE  B  OF 

FORM  APPROVED 

PAGES 

OMB  22  R  024S 

H.  EQUIPMENT  USAGE 


STAGE  1  _ _  STAGE  3 

STAGE  2 _ _  STAGE  4 


B.  TYPE  OF  SERVICE 


C  EXTENT  OF  USE 


D.  GEOGRAPHICAL  AREA  FOR 


STAGE  2 
STAGE  3 
STAGE  4 


E  NUMBER  OF  EQUIPMENTS  IN  INITIAL  STAGE 

F.  NUMBER  OF  EQUIPMENTS  PLANNED  FOR  OPERATIONAL  USE 

G  NUMBER  OF  EQUIPMENTS  OPERATING  SIMULTANEOUSLY  IN  THE  SAME  ENVIRONMENT  ^ 

H  TRANSMITTER  INSTALLATION 

1.  RECEIVER  INSTALLATION 

J  PREVIOUS  J/F  12  APPLICATION  NUMBER  TO  BE 

CONTINUED  UNCHANGED  J/F  12/ 

□  SUPERSEDED  J/F  12/ 

□  "*  LATE  D  J/f  >2/ 

|  |  NONE 


K 


L 


IS  AN  ECP  SPECTRUM  SIGNATURE  BEING  SCHEDULED  FOR  THIS  EQUIPMENT  IN  ACCORDANCE  WITH  MEASUREMENT  PROCEDURES 
DESCRIBED  IN  MiL-STD-449  I  )7  Q  y£s 


NO 


DESCRIBE  BRIEFLY  ANY  OTHER  MEASURED  DATA  BEING  OBTAINED  IN  LIEU  OF  OR  IN  ADDITION  TO  PROCEDURES  SET  FORTH 
IN  MIL -STD-449  (  I 


M  EQUIPMENT  SPECIFICATION  NUMBERISI 


iN  APPLICABLE  SERVICE  CONTRACT  NUMBERISI 


O  NAME  ANO  OFFICE  SYMBOL  OF 

PROGRAM  MANAGER 
PROJECT  ENGINEER 


TELEPHONE  NUMBERS 


Comrr>«rC'«l  _ 
Commtfcill  „ 


P  REMARKS 


CLASSIFICATION 


2.4  SIGNAL  DESIGN  AND  PROCESSING 


2.4.1  Introduction 

An  objective  of  the  Phase  I  program  has  been  to  consider  means  for 
optimizing  diversity  reception  of  VHF  communications.  Signal  processing 
can  provide  a  means  for  substantially  Improving  percent  copy  In  diversity 
communications  systems  by  producing  a  reliable  signal  from  multiple  signals 
that  may  be  unusable  Individually.  In  Phase  I  we  have  devoted  considerable 
attention  to  this  matter,  which  promises  to  be  of  significant  utility  In 
applications  of  diversity  reception  of  non-cooperative  VHF  transmitters. 
Although  the  Phase  II  program  emphasizes  the  resolution  of  outstanding 
physics  Issues  rather  than  signal  processing  Issues,  we  feel  that  It  Is 
Important  to  summarize  our  work  In  signal  design  and  processing  here,  so 
that  It  might  be  available  for  future  development  of  diversity  reception 
communications  systems. 

It  Is  well  known  that  the  Ionization  trails  produced  by  small  meteors 
as  they  disintegrate  In  the  high  atmosphere  result  In  short-term  communi¬ 
cation  channels,  which  are  usable  In  the  spectral  region  from,  roughly,  30 
to  100  MHz.  Since  the  formation  of  these  trails  Is  a  random  process,  the 
characteristics  of  the  meteor  burst  communication  link  between  two  given 
points  are  also  random  processes.  As  a  simplification,  we  may  think  of  the 
channel  as  being  either  on  or  off  with  typical  behavior  being  a  few  hundred 
"on”  periods  per  hour,  with  each  period  being  a  few  hundred  milliseconds, 
and  with  a  representative  duty  cycle  being  a  few  percent.  The  actual  value 
depends  on  the  characteristics  of  the  sender  and  the  receiver,  the  time  of 
year  and  of  day,  geographic  position,  etc. 

Whl  le  these  duty  cycles  permit  useful  traffic  transfer  between  co¬ 
operating  stations,  any  workable  system  for  monitoring  must  have  only  a  few 
percent  of  "off"  periods.  Such  a  large  Increase  In  duty  cycle  may  be 
achieved  by  forming  parallel  observation  paths  and  combining  signals  from 
these  multiple  paths.  That  Is,  several  antenna/receivers  may  be  oriented 
towards  the  location  to  be  monitored  and  their  output  relayed  to  a  central 
signal  processing  location  which  analyzes  them  Jointly  to  produce  an 
optimum  estimate  of  the  content  of  the  monitored  signal. 


Two  methods  of  providing  multiple  observation  paths  are  conceivable: 

e  Angular  diversity  where  many  co-located  narrow  beam  antennas  (plus 
a  receiver  per  antenna)  partition  the  volume  of  meteor  occurrences 
Into  several  high  gain  regions. 

•  Spatial  diversity  where  several  distributed  stations  with  moderate- 
beamwldth  antennas  provide  signals  from  distinct  volumes  of  space. 

The  two  approaches  are  Illustrated  In  Figure  2.22.  The  concept  of 
spatial  diversity  Is  self  evident.  Since  angular  diversity  partitions  the 
regions  of  meteor  trails*  each  beam  sees  fewer  events  of  a  given  strength. 
However*  the  Increased  gain  In  each  beam  permits  lower  Intensity  meteor 
trails  to  be  distinguished  (from  the  noise  floor)  so  that  the  duty  cycle  of 
the  narrower  beams  Is  affected  by  conflicting  factors.  We  assume  here  that 
It  remains  constant  as  the  beam  Is  narrowed. 

In  either  approach  then*  we  have  available  multiple  paths  that  are  at 
least  partially  Independent  In  terms  of  the  signal  propagation  mechanism 
and  are  essentially  Independent  In  terms  of  noise  Impairment.  The  signals 
are  Independent  when  they  arise  from  distinct  meteor  trails  and  are 
dependent  when  a  single  trail  gives  rise  to  a  receivable  signal  In  multiple 
paths.  Both  phenomena  give  rise  to  Increased  throughput.  The  Improvement 
Is  self  evident  for  Independent  signals.  With  dependent  signals  the 
additional  signal  power  from  multiple  observations  can  be  exploited  to 
extend  the  usable  duration  of  the  burst.  These  Improvements  require  that 
the  receivers'  direct  observables  be  brought  to  a  central  processing 
facility  so  that  the  combining  can  be  optimized.  In  the  angular  diversity 
case  such  joint  processing  Is  readily  achieved  since  all  receivers  are  co¬ 
located.  With  spatial  diversity*  however*  some  form  of  a  collection  net¬ 
work  Is  required  to  bring  the  various  observables  together  for  joint 
processing. 

The  signal  design  and  processing  methods  described  In  the  following 
are  viable  with  either  diversity  method.  The  combining  procedure  will 
allow  not  only  the  detection  of  all  signals  that  could  have  been  detected 
Individually*  but  also  may  permit  the  combination  of  Individually  unusable 
signal  remnants  Into  a  composite  signal  that  can  be  used. 

2.4.2  Channel  Description 

For  applications  of  a  diversity  reception  communication  system* 
multiple  propagation  paths  may  be  used  to  monitor  a  single  sender.  Each 
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Figure  2.22.  Diversity  Strategies  for  VHF  Monitoring 
via  Meteor  Bursts 


such  path  Includes  the  transmitter  being  monitored,  a  meteor  burst 
transmission  channel  and  a  receiver  (Including  an  antenna).  A  diagram  of 
the  overall  monitoring  network  for  certain  applications,  such  as  non- 
cooperative  reception.  Is  shown  In  Figure  2.23. 

These  multiple  paths  will  each  yield  an  observation  process  of  the 
sender’s  communication  activity  which  has  a  low  to  moderate  duty  cycle.  If 
the  paths  are  configured  to  yield  more  or  less  Independent  responses,  then 
the  joint  output  of  the  ensemble  of  receivers  may  be  processed  to  give  a 
diversity  combined  output  which  includes  the  combined  parallel  propagation 
effects  of  the  multiple  paths. 

If  these  paths  are  Independent,  then  each  path's  propagation  Is  the 
result  of  distinct  meteor  events.  Returning  to  the  "on-off”  conceptual 
model  of  the  last  section.  It  is  evident  that  a  parallel  path  monitoring 
network  would  result  In  a  combined  channel  which  Is  the  logical  OR  of  the 
set  of  Individual  channels.  Let  the  probability  that  any  given  path  Is  in 
the  "on"  state  be  PI.  Then,  assuming  independent  paths  and  identical 
Individual  duty  cycles,  the  probability  that  the  set  of  N  combined  paths  Is 
In  the  "on"  state  Is  PN  =  1— ( 1— P 1 )N .  Expanding  this  with  the  binomial 
theorem  yields: 

PN=  I^Pl" 

The  union  bound  of  probabil  Ity  theory,  PN  <  N*P1»  suggests  that  the 
duty  cycle  of  the  combined  channel  is  no  better  than  N  times  the  duty  cycle 
of  any  one  path.  Moreover,  even  this  level  of  Improvement  depends  strongly 
on  the  channels  being  statistically  Independent.  In  the  extreme  case  where 
the  same  meteor  phenomenon  Is  the  causal  mechanism  for  every  path,  then  the 
performance  of  the  combined  parallel  path  Is  no  better  than  the  performance 
of  any  one  path. 

The  results  of  the  previous  paragraph  are  unduly  pessimistic,  however, 
since  the  behavior  of  each  meteor  burst  channel  Is  more  complex  than  simple 
"on-off"  conditions.  Noise  Is  present  at  the  Input  to  each  of  the 
receivers  In  Figure  2.23  so  that  each  channel  would  be  "on"  only  In  a 
probabilistic  sense,  even  when  a  meteor  burst  exists.  A  more  precise  view 
Is  that  each  receiver  produces  a  continuous  random  variable  "r(t)".  If  we 
Ignore  for  the  moment  the  effects  of  variations  In  the  length  of  the 
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propagation  path*  then  r(t)  will  be  a  scalar  variable  which  Is  the  sum  of  a 
noise  process  and  of  a  signal  whose  amplitude  Is  a  random  process  generated 
by  meteor  bursts. 

We  assume  In  the  following  that  the  receiver  Input  noise  has  a  uniform 
spectrum  over  the  frequency  band  of  Interest  and  that  It  has  a  Gaussian 
distribution  so  that  an  additive  white  Gaussian  noise  (AWGN)  description  Is 
appropriate.  The  observation  process  from  each  of  the  receivers  Is  then  of 
the  form: 

r(t)=s(t)+w(t) 

where  "w(t)M  Is  white  Gaussian  noise  with  a  mean  of  zero  and  a  (one-sided) 
power  spectral  density  of  N0.  The  signal  process  ”s(t)M  consists  of  the 
transmitter’s  baseband  signal  "(ntt)”  modulating  a  carrier  of  angular  fre¬ 
quency  as  seen  after  passage  through  a  meteor  burst  channel.  This 
channel  may  be  characterized  by  Its  Instantaneous  gain  "git)"  and 
propagation  delay  ”d(t)”.  Thus  the  signal  observed  by  a  receiver  has  the 
form: 

s(t)  =  g(t)m[t-d(t)]cos[oo (t-d(t) )] 

At  first  glance  the  gain  process*  g(t)»  exhibits  the  typical  meteor 
burst  behavior.  That  Is*  It  Is  zero  much  of  the  time  (l.e.  burst  absent) 
and  exhibits  random  (Poisson  distributed)  bursts  of  finite  gain.  These 
bursts  usually  have  a  very  sharp  onset*  corresponding  to  the  rapid 
formation  of  an  Ionization  trail*  followed  by  an  exponential  decay  of  gain 
as  the  cylinder  of  Ionization  diffuses.  This  rate  of  decay  Is  reasonably 
small  compared  to  the  symbol  rate  of  transmission  of  the  source  being 
monitored  so  that  the  function  g(t)  may  be  approximated  as  a  constant  over 
the  symbol  Interval  for  symbols  of  the  order  of  10  milliseconds  or  shorter. 
The  delay  process,  d(t),  may  be  considered  a  constant,  to  a  first 
approximation,  during  the  course  of  a  given  meteor  burst  (l.e.  a  non-zero 
g(t)  episode).  The  delay  will  vary  from  one  burst  to  another*  however, 
depending  on  the  actual  location  of  the  Ionization  trail.  Consider  for 
example*  a  receiver  which  Is  1000  km  from  the  transmitter.  Then  a  meteor 
burst  occurring  at  an  altitude  of  100  km  above  the  mid-point  between  the 
transmitter  and  the  receiver  results  In  a  propagation  delay  of  about  3.40 
ms.  For  a  burst  at  the  same  altitude  but  above  a  point  200  km  at  right 
angles  from  the  mid-point*  the  transit  time  increases  to  3.65  or  about  0.25 
milliseconds  more.  For  symbols  of  the  order  of  10  ms  duration*  the 
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propagation  time  may  therefore  be  treated  as  constant  without  serious  error 
In  symbol  timing  (error  of  the  order  of  a  few  percent).  Propagation  time 
variations  of  a  few  hundred  microseconds  are  very  significant  on  the  time 
scale  of  the  carrier  signal#  however#  since  this  time  Interval  spans  a  few 
thousand  cycles  of  the  carrier.  From  another  perspective#  the  wavelength 
of  a  carrier  In  the  30-100  MHz  range  Is  a  few  meters  while  the  bursts  occur 
over  a  region  that  Is  hundreds  of  kilometers  In  extent  so  the  path  length 
of  two  different  meteor  bursts  may  be  expected  to  exhibit  a  carrier  phase 
which  varies  randomly  over  360  degrees.  There  may  also  be  some  variation 
In  the  path  length  (and  hence  of  the  carrier  phase)  during  the  course  of 
the  burst  as  a  result  of  displacement  of  the  Ionization  trail  by 
Ionospheric  winds.  Such  motion  results  In  the  received  carrier  appearing 
to  have  a  slightly  different  frequency  than  the  sender's  -  for  wind  speeds 
of  less  than  70  mph  (30  meter/second)  the  proportional  frequency  shift  will 
be  0.1  parts  per  million  or#  equivalently#  a  few  Hertz.  This  degree  of 
frequency  uncertainty  In  the  received  burst  signals  can  be  removed  by 
carrier  tracking  of  the  Incoming  signal  phase  so  long  as  the  symbol  period 
Is  short  enough  that  the  phase  change  is  small.  (This  will  be  the  case  for 
symbol  rates  of  a  few  tens  of  Hertz  or  higher). 

Putting  all  this  together  we  arrive  at  a  first  view  of  the  received 
signal  as: 

s(t)  =  gnm(t)cos(oot+  0n )  nT  <  t  <  (n+l)T 

The  above  expression  denotes  the  signal  during  the  nth  symbol  Interval  and 
states  that  the  channel  gain  Is  a  random  variable  which  Is  treated  as 
constant  over  the  symbol  Interval  and  that  the  signal  Is  received  with  a 
random  carrier  phase.  For  simplicity#  the  gross  delay  Is  Ignored  and  the 
expression  Is  written  In  the  time  frame  of  the  receiver. 

This  model  of  the  signal  process  Is  simplistic  In  the  sense  that  It  Is 
premised  on  either  a  single  relatively  strong  burst  or  no  burst  at  all. 
While  these  strong  bursts  constitute  the  usual  dominant  transfer  mode  for  a 
meteor  burst  channel#  the  nature  of  the  meteor  phenomena  Is  such  that  there 
Is  also  a  background  of  weaker  but  much  more  frequent  meteor  trails.  A 
more  complete  model  of  the  network  which  accounts  for  this  additional 
background  signal  activity  would  then  depict  a  scattered  background 
activity  with  a  Rayleigh  amp!  Itude  and  a  fluctuating  phase.  This  back¬ 
ground  Is  Interspersed  with  occasional  episodes  of  much  stronger  amplitude 


and  coherent  phase.  In  the  context  of  the  network  model  of  Figure  2.23# 
the  Incoherent  background  may  be  of  particular  Importance  since  the  back¬ 
ground  signals  In  the  multiple  channels  may  be  comblnable  to  yield  a  useful 
monitoring  of  the  transmitter  during  Intervals  when  no  major  burst  Is  In 
progress. 

It  should  be  emphasized  that  the  proceeding  signal  model  Ignores 
variations  In  the  transmltter-to-recelver  transit  time.  As  argued  above, 
this  Is  justifiable  In  the  context  of  low  symbol  rate  transmission.  How¬ 
ever,  for  monitoring  of  a  higher  speed  data  transmitter  (e.g.  4800  bps  or 
more),  the  time-dispersion  resulting  from  path  length  variation  will 
present  additional  signal  processing  chal  1  enges  to  deskew  the  multiple 
transmission  paths.  The  signal  processing  methodology  described  in  the 
fol  lowing  avoids  this  compl  icating  factor  by  addressing  only  low  speed 
symbols.  This  limitation  of  the  monitoring  method  should  be  kept  In  mind. 

In  earlier  paragraphs,  the  issue  of  channel  Independence  was  raised  In 
the  context  of  the  "on-off”  channel  model.  In  that  simple  context,  signals 
in  multiple  channels  from  a  common  meteor  could  not  enhance  the  duty  cycle 
above  that  for  the  single  channel.  Commonly  observed  signals  can  however, 
be  processed  to  yield  a  superior  duty  cycle  In  the  more  complete  model 
since  (largely)  Independent  noise  -  If  not  signal  -  Is  present  at  each 
receiver’s  Input  and  since,  with  proper  processing,  the  signals  from  the 
parallel  paths  may  be  combined  coherent! y.  This  results  In  a  stronger 
signal  which  will  have  a  longer  lifetime  before  dropping  beneath  the  noise 
floor. 

The  performance  of  the  network  model  may  then  be  summarized  as 
follows: 

•  The  multiplicity  of  paths  provides  for  the  observation  of  more 

meteor  burst  trails  since  a  larger  region  of  space  Is  observed 
(spatial  diversity)  or  several  regions  are  observed  with  enhanced 
gain  (angular  diversity). 

•  The  multiplicity  of  paths  may  permit  the  Incoherent  background 
activity  of  low  level  meteorites  to  be  exploited  for  enhanced  duty 
cycle. 

•  The  multlpl  Iclty  of  paths  can  result  In  an  array  gain  for  those 
meteor  trails  which  are  observed  In  more  than  one  beam. 


2.4.3  Experimental  Signal  Set 

From  the  discussion  of  the  channel  model  above*  It  Is  evident  that  the 
experimental  signal  set  should  have  two  key  properties.  First  It  should 
have  a  low  symbol  rate  to  avoid  the  time  dispersion  difficulties  alluded  to 
In  the  above  section.  Second*  It  should  be  robust  with  respect  to  the 
random  phase  characteristic  of  the  channel. 

Binary  frequency  shift  keying  at  100  symbols/second  Is  the  recommended 
signal  set  for  an  experiment  that  would  address  signal  processing  Issues. 
The  symbol  Interval  of  10  mS  Is  long  with  respect  to  the  dispersion  of 
propagation  time  so  this  signal  set  satisfies  the  first  criterion.  If  the 
two  frequencies  of  a  binary  signal  pair  are  each  Integer  multiples  of  the 
symbol  rate*  then  the  two  signals  In  the  set  are  orthogonal  over  the  symbol 
time.  An  orthogonal  signal  pair  may  be  distinguished  Incoherently  (l.e. 
without  knowledge  of  the  carrier  phase)  so  this  signal  set  also  satisfies 
the  second  criterion.  Such  a  set  can  also  be  processed  coherently*  of 
course,  with  significant  performance  advantage  under  low  contrast 
conditions. 

We  propose  then,  that  the  signal  set  be  comprised  of  100  and  200  Hz 
tones  each  of  duration  10  ms.  Representative  waveforms  In  this  class  are 
shown  In  Figure  2.24.  Mathematically*  this  signal  set  Is  defined  as: 
mn(t»j)  *  As1n(27rfj  (t-nT)) 

where: 

n  denotes  the  symbol  epoch  (n=0»l*2*  ...) 

T  Is  a  symbol  Interval  (10  milliseconds) 

Jn  represents  the  choice  of  frequency  for  epoch  n  and  Is  binary  (0 
or  1). 

f  Is  the  Instantaneous  frequency  during  the  symbol  Interval*  and  Is 
one  of  two  choices  fQ=100  Hz  or  f^=200  Hz. 

A  Is  the  amplitude  of  the  signal  set  and  Is  a  constant. 

Fortuitously,  a  signal  set  with  these  properties  seems  to  be 
representative  of  operational  traffic  of  interest  for  monitoring  via  a  VHF 
reception  network. 

The  proposed  signal  set  has  certain  other  advantages  for  experiment. 
The  low  symbol  rate  results  In  a  relatively  large  energy-to-nolse  ratio  for 
moderate  carrier  powers  so  that  many  meteors  should  give  rise  to  usable 
signals.  Concommltantly*  the  low  rate  signal  has  a  narrow*  concentrated. 
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spectrum  which  minimizes  out-of-band  spectral  emissions.  The  constant 
amplitude  waveform  of  these  signals  will  minimize  the  consequences  of  any 
nonl Inearltles  In  amplifiers  through  which  the  signal  Is  processed.  Of 
considerable  Importance  to  the  experiment*  these  low  rate  signals  lend 
themselves  to  processing  by  general  purpose  micro-computers*  thus 
minimizing  program  cost  and  maximizing  program  flexibility. 

The  data  content  of  the  signal  depends  on  the  stream  of  binary  letters 

(the  ”jn"  fn  the  equation  above).  For  the  experiment  this  data  stream 
should  be  a  clearly  definable  pattern  so  that  the  outputs  of  Individual 
receivers  and  of  the  central  diversity  combining  can  be  readily  checked  for 
correctness.  On  the  other  hand*  the  transmitted  data  pattern  should  appear 
more  or  less  random  so  that  good  balance  occurs  between  logical  0  and  1  and 
so  that  the  potential  for  biases  In  the  signal  processing  Is  minimized. 
The  so  called  psuedo-nolse  (PN)  sequences  have  these  properties  and  are 
very  easy  to  Implement.  Recommended  for  a  signal -processing  experiment  Is 
a  PN  sequence  of  periodicity  8,388,607  which  may  be  generated  by  a  linear 
feedback  shift  register  of  length  23.  A  suitable  recursion  relation  to 
generate  a  PN  sequence  of  this  type  Is: 

Sn+23  =  Sn+Sn+3  (modulo  2  addition) 

Note  that  the  periodicity  of  this  sequence  Is  very  nearly  one  day  for  a  bit 
rate  of  100  per  second. 

2.4.4  Modulator  Organization 

The  modulator  will  generate  the  PN  data  sequence  and  map  this  sequence 
Into  the  frequency  shift  keyed  waveform  described  above.  These  operations 
will  be  Implemented  In  a  smal  1  general  purpose  computer  equipped  with  a 
digital -to-analog  convertor.  The  technique  will  be  to  produce  a  quantized 
version  of  the  FSK  waveform  using  a  sampled  construction  which  Includes 
many  samples  per  symbol  of  the  experimental  waveform.  A  sample  rate  of  16 
times  or  more  the  bit  rate  Is  appropriate  and  easily  realizable.  After 
passage  through  low  pass  filtering  to  remove  the  high  frequency  artifacts 
of  the  sampled  reconstruction,  the  resultant  signal  will  closely  follow  the 
signal  set  described  above.  Note  that  this  signal  set  Is  equivalent  to  a 
single  sinusoidal  cycle  per  symbol  or  two  such  cycles  and  that  exact  lock 
will  be  maintained  between  the  symbol  timing  and  the  FSK  tone  frequencies 
so  that  al 1  transitions  are  continuous. 


The  modulator's  output  will  be  presented  to  the  VHF  transmitter  which 
operates  single  sideband.  Thus,  the  actual  transmitted  waveform  will  also 
be  In  constant  envelope  FSK  format  with  100  Hz  separation  between  Its  tones 
and  with  a  100  Hz  symbol  rate. 

A  top  level  block  diagram  of  the  modulator  hardware  In  shown  at  Figure 
2.25.  A  WWV  Interface  has  been  posited  In  the  diagram  as  a  means  of 
locking  the  modulator  to  a  standard  timing  source.  The  structure  of  soft¬ 
ware  to  generate  the  experimental  signals  Is  shown  In  Figure  2.26.  Note 
that  In  this  and  other  software  diagrams  to  follow.  It  has  been  assumed 
that  the  transmitter  and  receivers  operate  under  the  control  of  a  central 
processing  site  via  telephone  lines.  Some  form  of  packet  protocol  has  been 
assumed  as  the  means  of  exchanging  Information  on  this  wirel ine  system. 
Provisions  are  also  Indicated  In  Figure  2.26  for  a  man-machine  Interface  at 
the  transmitter  equipment. 

2.4.5  Demodulator  Organization 

The  use  of  a  microcomputer  for  signal  processing  Is  a  feature  of  the 
demodulator  also.  A  signal  processing  Interface  consisting  of  A/D  and  D/A 
convertors  Is  an  essential  adjunct  to  the  microcomputer  for  this  purpose. 

Before  exploring  the  actual  Implementation  structure,  however,  it  Is 
worthwhile  here  to  further  consider  some  fundamental  signal  processing 
Issues.  When  a  meteor  burst  occurs,  the  signal  arriving  at  the  demodulator 
has  a  random  phase.  While  It  Is  not  essential  to  recover  the  carrier  phase 
to  demodulate  an  FSK  signal,  there  Is  a  performance  advantage  to 
demodulating  coherently  where  possible.  We  propose  to  preserve  the 
potential  for  coherent  detection  In  the  demodulator  design  by  projecting 
the  received  signal  onto  quadrature  components  (I  &  Q)  as  the  first  proces¬ 
sing  step  In  the  demodulator.  Since  this  Is  a  reversible  step,  no  loss  of 
optimality  results  even  If  the  exact  phase  relation  Is  unknown.  (There  Is 
an  assumption,  however,  that  the  carrier  frequency  of  the  transmitter  Is 
known  within  a  small  uncertainty  so  that  the  projection  of  the  signal  onto 
the  phase  space  of  the  receiver  may  be  done  without  rapid  change  of  phase). 

With  reference  to  Figure  2.27,  which  shows  a  top  level  block  diagram 
of  the  demodulator,  the  resolution  of  the  I  and  Q  components  Is  performed 
In  the  block  labelled  "audio  synchronous  mixers".  The  structure  of  this 
block  appears  In  more  detail  In  Figure  2.28,  which  shows  that  the  operation 
consists  of  two  down-conversion  frequency  translations  using  orthogonal 


Figure  2.25.  Modulator  Hardware  Organization 


received  observables  on  a  long-term  storage  medium  (magnetic  tape)  Is  made 
so  that  the  data  may  be  processed  using  advanced  combining  algorithms  If 
and  when  they  become  available. 


obtained.  Such  Intermediate  level  signals  may  still  be  detected 
Incoherently  on  the  basis  of  the  relative  energy  in  the  correlators  for 
each  of  the  two  tones.  (That  is*  by  forming  the  sum  of  the  square  of  the  I 
and  the  Q  correlators  for  each  frequency).  The  CPS  may  then  combine  the 
energy  estimates  for  each  frequency  from  each  of  the  monitoring  network's 
parallel  paths  into  a  combined  decision  variable  for  each  of  the  two  tones. 
This  diversity  combining  should  weight  the  various  components  according  to 
the  current  quality  of  meteor  propagation  for  the  various  paths.  Thus,  the 
observable  from  a  stronger  path  would  receive  more  credance  than  that  from 
a  weaker  path.  Such  weighting  requires  a  procedure  for  estimating  the 
relative  quality  of  moderate  to  weak  meteor  burst  propagation  paths, 
however.  Both  this  channel  assessment  procedure  and  the  combining 
algorithm  require  further  system  analysis  as  well  as  software  development 
and  hence  are  appropriate  for  some  later  program  exploring  signal 
processing  Issues. 

The  top  level  organizations  of  the  CPS  hardware  and  software  are  shown 
In  Figures  2.32  and  2.33,  respectively.  Packets  containing  the  various 
paths'  received  observables  will  arrive  at  the  CPS  from  the  telephone 
network  or  from  direct  connection  with  the  network's  receivers.  Since  each 
of  the  propagation  paths  will  have  a  specific  average  delay,  the  components 
from  the  various  receivers  must  be  aligned  (or  de-skewed)  before  they  can 
be  combined.  (This  deskewing  should  not  be  thought  of  as  taking  out  the 
variation  In  any  one  path  due  to  differences  In  that  path's  length  from  one 
burst  to  the  next.  It  simply  removes  the  average  difference  between  the 
propagation  time  for  two  distinct  paths).  Also,  the  apparent  signal 
strength  must  be  assessed  so  that  the  appropriate  diversity  combining  and 
detection  algorithm  may  be  employed.  In  the  case  of  strong  signals,  the 
phase  trajectory  of  the  signal  Is  determined  beforehand  to  permit  coherent 
detection  to  be  done.  Detection  proceeds  by  whatever  algorithm  Is  most 
appropriate  to  the  quality  of  the  data  and  yields  a  composite  network 
decision  as  to  the  monitored  source's  data  bit.  The  relative  quality  of 
these  decisions  Is  appended  to  them  and  the  result  is  output  to  a  network 
operator. 

The  CPS  will  be  fitted  with  additional  software  to  permit  an  operator 
to  control  remote  recel vers.  Additional  features  permit  an  operator  to 
review  the  overall  performance  of  a  network.  Provision  to  record  all  the 


composite  observable  (l.e.  six  vectors  of  four  components  each).  In  the 
terminal ogy  of  detection  theory#  a  maximum  aposterlor  probabl  1  Ity  (MAP) 
decision  Is  required.  For  equally  likely  0’s  and  l’s  at  the  transmitter, 
this  requirement  readily  maps  (through  Bayes’  Rule)  to  an  equivalent 
requirement  for  a  maximum  likelihood  (ML)  decision.  That  Is#  the  CPS’s 
task  Is  to  find  the  binary  value  which  satisfies 

Max  Pt^l#£2#r3»I14»£S*£6  !  Jn3 

jn  =  0  or  1 

With  reference  to  the  subsection  on  demodulation#  above#  It  may  be  recalled 
that  the  four  components  of  each  of  these  vectors  are  comprised  of  a  signal 
component  plus  a  noise  component.  The  noise  may  be  treated  as  Independent 
and  Identically  distributed  Gaussian  random  variables.  The  signal 
components  of  the  two  correlators  which  are  matched  to  the  tone  that  was 
actually  sent  have  means  of  the  form  gnAcos(0n)  and  gnAs1n(0n)  while  the 
other  two  components  have  means  which  are  zero.  Both  gn  and0n  are  random 
variables  (as  a  result  of  meteor  burst  activity)  but  may  be  treated  as 
constants  over  the  period  of  one  symbol. 

From  a  practical  point  of  view  It  may  be  difficult  to  make  decisions 
In  a  space  of  this  high  dimensionality.  A  sub-ideal  but  more  readily 
Imp  1 ementab  1  e  approach  Involves  pre-processing  each  receiver’s  signal 
vector  to  establish  If  any  of  them  are  experiencing  a  strong  meteor  burst. 
When  such  a  burst  Is  In  progress  on  a  path#  this  path  will  be  the 
predominate  factor  In  reaching  a  network  decision. 

To  keep  Initial  software  development  from  becoming  too  complex#  the 
Initial  phase  of  the  CPS  should  use  this  approach  to  signal  processing. 
That  Is#  It  should  assess  which  of  the  six  receiver  inputs  have  strong 
bursts  In  progress  and  Ignore  the  rest.  For  those  with  strong  bursts,  the 
CPS  should  estimate  the  carrier  phase  and  phase  rate  for  the  current 
propagation  path.  With  a  reliable  phase  estimate,  the  sequence  of  four 
dimensional  observations  from  the  receiver  can  be  reduced  to  two 
dimensions.  If  more  than  one  channel  Is  responding  to  a  strong  burst  In  a 
given  symbol  Interval  then  these  two  dimensional  representations  may  be 
added  coherently  to  produce  a  better  equivalent  signal -to-nolse  ratio  prior 
to  detection. 

At  the  next  level  of  complexity  the  CPS  should  be  enhanced  to  process 
weaker  signals  which  are  not  sufficiently  strong  for  a  phase  estimate  to  be 


of  angular  diversity*  the  receivers  may  all  be  situated  In  the  same  loca¬ 
tion  so  this  Interconnection  may  be  very  straightforward  -  a  direct  RS232 
Interface  for  example. 

For  spatial  diversity*  Individual  receivers  may  be  scattered  over  a 
wide  region.  An  approach  to  Interconnection  In  this  case  Is  to  use  the 
public  telephone  network  with  medium  speed  wireline  modems  at  the  receiver 
sites  and  the  central  processing  site.  The  observable  vectors  are 
transferred  from  the  receivers  to  the  central  processor  as  they  are 
received.  Recalling  that  the  signalling  rate  Is  100  per  second  and  that 
each  vector  contains  four  components,  there  are  400  such  components  per 
second  to  transfer  for  each  receiver.  If  the  sampling/processing  resolu¬ 
tion  of  16  bits  per  sample  Indicated  In  Figure  2.30  Is  preserved  for  the 
Interconnection,  then  the  resultant  data  rate  Is  some  6400  bits  per  second. 
This  data  transfer  rate  Is  above  the  limits  of  low  cost  wireline  modems  and 
of  unconditioned  voice  lines.  With  some  reduction  In  resolution  (e.g.  12 
bits/sample)  and  a  50%  discard  rate  by  the  weak  signal  detection  algorithm, 
the  rate  drops  to  2400  bits/second  which  may  be  handled  with  moderate  price 
equipment  and  lines. 

2.4.7  Central  Processor  Organization 

Once  all  the  signals  from  the  parallel  components  of  the  network  are 
brought  to  a  common  location  they  can  be  processed  jointly  to  achieve  some 
diversity  and/or  array  gain.  The  detailed  approach  to  Joint  processing 
will  require  further  development  during  the  next  phase  of  the  project. 
However,  some  general  Issues  are  evident  at  this  time  and  are  discussed 
here.  (For  the  sake  of  definiteness  in  the  discussion,  six  paths  are 
assumed.  This  number  may  be  greater  or  less  depending  on  the  detailed 
design  of  an  experiment.  Clearly,  at  least  two  paths  are  required). 

For  each  signalling  Interval  of  the  transmitter  being  monitored,  the 
central  processing  site  (CPS)  will  be  represented  with  six  observable 
vectors,  one  from  each  receiver,  with  each  vector  having  four  components. 
(See  the  preceding  sections  for  a  description  of  these  observables).  For 
observations  which  seemingly  have  no  useful  signal  content,  no  data  is 
forwarded  to  the  CPS  by  the  receiver  Involved:  the  CPS  will  treat  these  as 
pure  noise. 

In  an  Ideal  sense  the  problem  facing  the  CPS  Is  to  decide  which  of  the 
binary  hypothesis  Is  more  probable  to  have  been  transmitted  given  the 
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Figure  2.31.  Demodulator  -  Software  Organization 


as  by  oscillator  differences.  It  Is  Important  to  keep  the  corrective 
changes  to  the  local  oscillator  heavily  damped  to  avoid  unnecessary  Jitter 
In  the  frequency  tracking. 

The  other  major  housekeeping  function  Is  symbol  timing  tracking.  Note 
that  the  proposed  algorithm  for  maintaining  this  timing  depends  heavily  on 
the  relative  stability  of  the  transmitter's  and  the  receiver’s  symbol 
clocks  and  the  presence  of  an  occasional  strong  meteor  burst  to  refresh  the 
alignment  of  these  clocks.  To  assess  the  reasonableness  of  this  approach 
consider  that  It  Is  not  difficult  to  provide  a  clock  stability  of  1  ppm  at 
each  end  of  the  path.  At  a  symbol  rate  of  100  symbols/second  (360,000  per 
hour)  the  maximum  slip  rate  Is  less  than  one  symbol  position  per  hour. 
Since  there  are  typically  many  tens  to  a  few  hundred  strong  meteor  bursts 
per  hour,  these  should  provide  ample  opportunity  for  correction  of  the 
receiver's  symbol  timing  to  within  a  small  fraction  of  a  symbol  Interval. 
As  for  AFC,  the  symbol  timing  correction  algorithm  Is  engaged  only  during 
the  presence  of  strong  signals  and  Is  decision  directed.  The  algorithmic 
basis  will  be  an  early-1  ate  gate  which  compares  the  relative  content  of 
each  tone  (100  or  200  Hz)  In  a  window  which  the  algorithm  attempts  to 
center  on  the  transition  time  from  one  symbol  to  a  next. 

As  noted  earlier,  a  symbol  period  is  assumed  to  be  long  relative  to 
the  variation  in  the  propagation  time  as  a  result  of  differing  path 
lengths.  Such  Is  the  case  for  the  100  bit  per  second  signals  proposed  for 
the  experiment.  A  more  complex  symbol  timing  algorithm  will  be  required, 
however.  If  the  signalling  speed  Is  made  substantially  higher. 

The  organization  of  software  to  perform  the  demodulator  tasks  within  a 
microcomputer  Is  shown  In  Figure  2.31.  Many  of  the  blocks  of  this  diagram 
relate  directly  to  the  preceding  discussion  of  the  demodulator.  Additional 
blocks  Include  provision  for  organizing  the  observable  vectors  Into  packets 
for  transmission  to  a  central  processing  site,  for  man-machine  Interac¬ 
tions,  for  recording  the  observed  data  on  mass  storage  (e.g.  a  hard  disc), 
and  for  keeping  the  tlme-of-day  In  unison  with  the  central  processing  site. 
2.4.6  Interconnection  of  Receivers 

For  either  of  the  two  types  of  diversity  discussed  In  this  subsection, 
there  Is  a  need  to  bring  together  the  receiver  output  signals  from  the 
multiple  paths  to  a  common  location  for  diversity  combining.  In  the  case 
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It  may  be  desirable  to  limit  the  transfer  of  receiver  observables  to  those 
which  appear  to  be  at  least  marginally  useful.  Thus  each  receiver  will  run 
a  weak-signal -present  algorithm  to  differentiate  between  unusable  noise 
outputs  and  those  with  sufficient  signal  remnant  to  warrant  transfer  to  the 
central  processing  site.  The  basis  of  the  algorithm  will  be  the  ratio  of 
the  apparent  relative  energy  In  the  100  Hz  correlators  to  that  in  the  200 
Hz  correlators. 

The  operation  of  the  demodulators  may  then  be  summarized  as  being  a 
set  of  four  correlators  which  provide  matched  filtering  for  each  of  the  two 
tones  In  each  of  the  quadrature  phases  of  the  received  carrier.  These  yield 
a  four  component  vector  observable  at  each  symbol  time.  Preliminary 
processing  of  the  observables  Is  done  to  cull  noise-only  outputs  from 
possible  signal  present  observables  which  are  forwarded  to  a  central 
processing  site.  Housekeeping  functions  essential  to  demodulator  operation 
Include  correction  of  the  demodulator's  symbol  clock  to  keep  It  synchronous 
with  the  transmitter's  and  correction  of  the  demodulator's  local  oscillator 
frequency  to  compensate  for  frequency  differences  between  the  transmitter 
and  the  receiver. 

Both  of  these  housekeeping  algorithms  are  engaged  only  when  a  strong 
signal  component  Is  manifestly  present  -  l.e.  when  a  strong  meteor  burst  is 
In  progress.  (Detection  of  the  presence  of  strong  signals  Is  a  simple 
extension  of  the  weak  signal  algorithm  previously  discussed).  The  auto¬ 
matic  frequency  control  (AFC)  algorithm  will  be  based  on  the  rate  of  change 
of  phase  of  the  detected  bit  stream.  That  Is.  the  four  component 
correlator  output  will  be  comprised  of  a  strong  signal  component  in  the 
output  of  either  the  100  or  200  Hz  correlators  and  of  noise  only  in  the 
other  pair  of  correlators.  Suppose,  for  example,  that  a  100  Hz  tone  Is 
sent  In  the  nth  symbol  Interval  and  that  a  strong  burst  Is  In  progress. 
Then  the  outputs  of  the  100  Hz  correlators  In  the  I  and  Q  channels  will  be 
of  the  form  gnAcos(  0n)  and  gnAs1n<  0n>*  respectively  (plus  some  noise). 
The  ratio  of  the  Q  to  I  correlator  outputs  produces  an  estimate  of 
tan(0n).  From  several  such  estimates  in  succession  an  estimate  of  the 
rate  of  change  of  the  phase  angle  per  symbol  interval  can  be  made  and  used 
to  generate  a  corrective  signal  to  the  local  oscillator  to  bring  the  rate 
of  change  to  zero.  Since  the  apparent  frequency  of  any  Incoming  burst  Is 
Influenced  by  Doppler  shifting  from  motion  of  the  Ionization  trail  as  well 
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WHERE  s(n,0)  IS  THE  QUANTIZED  VALUE  OF 
sin  (2  TTn/16) ,  SUITABLY  SCALED. 


Figure  2.30.  Demodulator  -  Implementation  of  Correlators 


Ai  =  tl0»n»Il»n»Q0»n*dl»n) 

A1 1  the  Information  that  Is  useful  In  making  a  decision  as  to  the  trans¬ 
mitter's  nth  binary  Input  CI.e.  jn>  Is  preserved  In  this  observation 
vector.  For  a  one-path  observation  of  the  monitored  transmitter*  a 
decision  on  jn  would  be  made  on  the  basis  of  this  observation  vector  alone. 
For  an  acutal  signal -processing  experiment*  many  concurrent  observations  of 
the  transmission  will  be  available.  Optimum  decisions  for  the  network  as  a 
whole  require  that  these  observations  be  processed  jointly.  Hence  In  an 
Ideal  sense  It  Is  not  appropriate  for  each  receiver  to  Infer  a  decision  on 

jn  from  tn  so  each  receiver  will  simply  pass  along  these  observation 
vectors  to  a  central  processing  site. 

Several  remarks  are  apropos  here.  First#  It  should  be  clearly  noted 
that  an  Implicit  assumption  was  made  that  the  receiver  Is  aware  of  the  time 
frame  of  the  transmitter  to  a  sufficient  degree  of  accuracy  to  carry  out 
the  correlation.  (The  correlation  Integration  proceeds  across  the  span  of 
a  symbol  time  and  must  be  synchronously  aligned).  This  assumption  will  be 
justified  presently.  Second,  the  Ideal  correlation  operation  shown  In 
Figure  2.29#  Is  difficult  to  Implement  exactly.  A  close  approximation  Is 
shown  In  Figure  2.30  which  Implements  the  correlator  digitally  In  discrete 
time  and  discrete  amplitude.  This  digital  correlator  super-samples  the 
analog  waveform  with  a  quantizer  (A  to  D  converter)  to  an  appropriate 
resolution  (12  bits  Is  achievable  at  low  cost)#  weights  these  by  the 
template  of  the  prototype  signal  and  then  carries  out  a  simple  digital 
Integration.  (A  super  sampling  rate  of  16  times  the  symbol  rate  Is  assumed 
here:  the  exact  value  Is  flexible.)  Third#  there  are  housekeeping 
functions  that  the  demodulator  must  perform  which  are  best  Implemented  on  a 
"decision  -  directed"  basis.  (These  Include  symbol  timing  tracking  and 
automatic  frequency  control).  Thus  under  some  conditions#  the  demodulator 

will  further  process  rn  to  obtain  a  best  estimate  of  Jn*  This  estimate  Is 
only  used  locally  for  these  housekeeping  functions*  however#  and  the  net¬ 
work  estimate  of  Jn  Is  based  on  the  set  of  £n  ^rom  receivers.  Fourth# 
while  the  Ideal  requires  that  all  receiver's  observables  should  enter  Into 
a  network-wide  decision  variable#  the  reality  Is  that  signals  below  some 
minimum  level  are  of  negligible  value  In  reaching  a  composite  decision. 
Since  the  transport  of  this  data  from  a  remote  receiver  to  a  central 
processing  site  requires  communication  capacity  (e.g.  on  a  wireline  modem) 
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(sine  and  cosine)  variants  of  a  local  oscil lator  output.  It  should  be 
noted  here  that  the  SSB  transmitter  and  receiver  are  assumed  to  be  tuned  to 
slightly  different  center  frequencies  -  a  difference  of  1800  Hz  has  been 
assumed  for  discussion.  The  output  of  the  receiver  Is  then*  nominally#  an 
FSK  signal  with  frequencies  of  1900  and  2000  Hz.  The  audio  synchronous 
mixers  down  convert  to  the  original  FSK  process  (of  100  &  200  Hz)  as  an 
Integral  part  of  the  quadrature  projection  activity.  Note  also  that  the 
local  oscillator  In  the  audio  synchronous  mixers  Is  electrically  tunable. 
This  feature  Is  exploited  subsequently  in  the  demodulation  to  complete  an 
automatic  frequency  control  loop  which  Is  used  to  compensate  for  any  fre¬ 
quency  differences  between  the  VHF  radio  at  the  transmitter  and  the 
receiver. 

After  low  pass  filtering  to  remove  the  high  order  mixing  frequencies, 
the  quadrature  components  are  presented  to  the  demodulator’s  microcomputer 
where  they  are  sampled  and  correlated.  The  signals  have  thus  been  reduced 
to  their  original  baseband  form#  except  that  they  are  rotated  In  phase 
space  by  an  arbitrary  phase  angle  and  are  corrupted  by  additive  noise. 
Mathematically,  these  Inputs  are: 

rl(t,n)  *  gnAs1n(27rfjn  (t-nT) )cos(  Q  n  +  wj(t)) 

rQ(t#n)  =  gnAs1n(2  7rfjn  (t-nT) )s in (  6  n  +  WQ(t)) 

Here  the  representation  Is  for  the  nth  symbol  Interval  of  T  (  =  10)  milli¬ 
seconds#  the  frequencies  fjn  are  either  100  or  200  Hz#  the  unknown  phase 
angle  Is  0n  for  the  nth  symbol  time*  gn  Is  the  response  of  the  meteor 
burst  channel  for  this  Interval*  and  Wi(t)  and  WQ^t)  are  Independent  white 
noise  processes. 

It  Is  well  known  that  each  of  these  time  waveforms  may  be  converted  to 
binary  decision  scalars  by  correlating  the  waveform  against  each  of  Its 
possible  causal  antecedents  (l.e.  the  100  Hz  or  200  Hz  tone).  (Observe 
that  the  effect  of  the  unknown  phase  Is  to  simply  scale  the  signal 
magnitude  by  a  trigonometric  function.)  While  this  transformation  (cor¬ 
relation)  Is  not  reversible*  only  Irrelevant  portions  of  the  noise  are  lost 
In  the  operation.  The  correlation  operation  Is  diagrammed  at  Figure  2.29 
for  the  ”1”  channel.  An  Identical  structure  Is  required  for  the  ”Q” 
channel.  At  the  end  of  each  symbol  Interval*  there  are  two  scalar  outputs 
from  each  of  the  I  and  Q  correlators  to  yield  a  four  component  observation 
vector  which  we  denote  for  the  nth  symbol  Interval  as: 
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Figure  2.28.  Demodulator  Audio  Synchronous  Mixer 


Figure  2.27.  Demodulator  Hardware  Organ 1 zat 


Figure  2.26.  Modulator  Software  Organization 


Figure  2.32.  Central  Processor  -  Hardware  Organization 
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Figure  2.33.  Central  Processor  -  Software  Organization 
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3.0  TECHNICAL  APPROACH 


As  mentioned  In  the  Introduction#  the  VHF  Reception  Program  -  Phase  II 
Is  Intended  to  be  the  simplest  that  resolves  the  outstanding  physics  Issues 
concerning  diversity  reception  of  meteor  path  propagation.  These  Issues 
are  summarized  by  the  following  questions. 

•  Gain  vs.  Beamwldth 

What  Is  the  gain  vs.  beamwldth  tradeoff  on  a  single  receiver 

antenna? 

•  Angular  Diversity 

Can  Independent  meteor  propagation  paths  be  obtained  at  a  single 
location  using  multiple  receiver  antennas? 

•  Spatial  Diversity 

What  spatial  separation  Is  needed  to  obtain  independent  paths 
from  a  common  volume? 

We  have  found  that  the  simplest  and  least  expensive  approach  to 
addressing  these  Issues  Is  one  suggested  to  us  by  Dr.  Fred  Quelle#  the 
Scientific  Officer  of  the  Phase  I  program.  This  approach  calls  for  a 
central  receiver  site  and  multiple  transmitter  sites.  The  central  receiver 
site  contains  three  Yagl  antenna  arrays  and  a  single  receiver.  It  should 
be  located  where  the  Phase  II  program  will  be  principally  performed.  The 
transmitter  sites  should  be  located  at  appropriate  locations  to  achieve  the 
necessary  spatial  and  angular  diversity.  The  transmitters  should  be 
unattended#  and  turned  on  and  off  by  simple  telephone  commands. 

The  simplest  approach  has  also  been  taken  for  signal  processing.  To 
resolve  the  physics  Issues  of  Phase  II#  the  data  required  Is  percent  copy 
of  each  transmitter.  Therefore#  the  simplest  signal  from  each  transmitter 
Is  continuous#  steady  transmission  at  a  single  frequency.  To  differentiate 
the  signals  from  each  Individual  transmitter  during  simultaneous  transmis¬ 
sion#  each  transmitter  would  operate  at  a  slightly  different  sideband  of 
the  same  carrier  frequency#  such  as  40  MHz.  Four  transmitters  might 
transmit  at  40  MHz  +  500  Hz,  40  MHz  +  800  Hz,  40  MHz  +  1100  Hz,  and  40  MHz 
+  1400  Hz.  At  the  receiving  site#  the  principal  carrier  frequency  would  be 
subtracted  from  the  Incoming  signals#  so  that  the  transmitter  signals  would 


produce  distinct  audible  tones  at  500  Hz,  800  Hz,  1100  Hz,  and  1400  Hz. 
The  data  would  then  be  collected  on  a  tape  recorder  for  subsequent 
analysis. 

It  Is  worth  reemphasizing  that  the  objective  of  the  Phase  II  program 
Is  not  to  simulate  applications,  such  as  reception  of  noncooperative 
tactical  VHF  communications,  but  rather  to  address  physics  Issues  that  must 
be  resol  ved  In  order  to  devel op  the  capabll  Ity  to  design  optimal  netted 
diversity  reception  communications  systems.  This  objective  justifies,  for 
example,  the  use  of  steady  tones  as  signals,  and  justifies  Invoking 
reciprocity  to  consider  a  system  with  multiple  transmitters  and  a  single 
receiver,  even  though  reciprocity  does  not  strictly  hold  In  meteor  burst 
communications,  as  explained  In  Subsection  2.2.  The  duty  cycle  differs 
when  a  receiver  and  transmitter  exchange  locations,  owing  to  asymmetries  In 
the  number  of  meteor  bursts  over  different  locations,  particularly  along 
east-west  paths.  The  JAYC0R  MBCS  code  accounts  for  these  asymmetries, 
however,  and  provides  a  valid  standard  against  which  to  compare  data  from 
any  receiver/transmitter  configuration.  Therefore,  the  proposed  configura¬ 
tion  Is  just  as  valid  for  addressing  the  physics  Issues  as  the  more  compex 
and  expensive  alternatives  that  have  been  considered. 

One  disadvantage  of  the  proposed  configuration  Is  that  multiple  trans¬ 
mitters  Imply  multiple  FCC  authorizations,  and  the  authorizations  must  be 
coordinated  In  frequency.  We  have  been  told,  however,  that  obtaining 
authorizations  for  government-sponsored  transmissions  Is  a  relatively 
straightforward  procedure  requiring  only  90  to  180  days  In  many  cases. 
Since  the  authorizations  must  be  coordinated  In  frequency  and  are  site- 
specific,  somewhat  more  time  should  be  allowed  for  obtaining  the  authoriza¬ 
tions,  and  ample  redundancy  should  be  provided  In  the  number  of  sites  for 
which  one  seeks  authorizations. 

The  transmitter  sites  should  be  provided  by  military  bases  and  govern¬ 
ment  property  at  no  cost  for  land  use  to  the  Phase  II  program.  Indications 
are  that  military  bases  will  provide  land  and  utilities,  but  will  want  the 
experimenters  to  provide  a  shed  or  trailer  for  their  equipment. 

We  have  received  preliminary  approval  to  locate  a  central  receiving 
site  at  Camp  El  lot,  which  Is  In  the  eastern  part  of  NAS  Miramar.  Because 
Camp  Eliot  Is  at  the  northeastern  outskirts  of  San  Diego,  and  the  receiving 
antennas  would  be  looking  north  and  east  only.  Interference  by  RF  noise 
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from  San  Diego  Is  expected  to  be  moderate.  Housing  for  the  experimental 
equipment  will  also  be  available  at  Camp  Eliot. 

This  completes  the  discussion  of  the  strategy  of  the  overall  approach 
to  the  Phase  II  program.  The  rationale  behind  the  specific  technical 
approach  to  each  of  the  physics  Issues  Is  summarized  in  the  remainder  of 
this  section. 

3.1  Gain  vs.  Beamwldth 

To  a  good  approximation,  the  power  in  the  main  lobe  of  an  antenna  Is 
Inversely  proportional  to  the  angular  width  of  the  main  lobe.  This 
relation  allows  some  flexibility  In  the  design  of  antennas  for  meteor-burst 
communications  systems.  To  date,  virtually  all  meteor-burst  communications 
systems  have  used  wlde-beam  antennas,  with  main  lobes  having  3  dB  points 
typically  separated  by  about  40°  or  50°. 

An  unresolved  physics  Issue  Is  whether  It  may  be  advantageous  to  focus 
a  receiving  antenna  to  a  smal  ler  beam  at  a  higher  gain.  Of  course,  the 
smaller  beam  subtends  a  smaller  field  of  view  and  encompasses  fewer  meteor 
bursts  per  hour.  But  the  higher  gain  of  the  antenna  compensates  to  some 
extent  by  making  more  reflected  signals  observable  per  unit  solid  angle. 
Whether  there  might  be  a  net  gain  of  percent  copy  by  focusing  a  receiving 
antenna,  and  what  the  optimal  gain  and  beamwldth  are  for  a  given 
transmitter/receiver  configuration  are  the  issues  to  be  addressed. 

The  Monte-Carlo  code,  MBCS,  which  was  developed  In  the  Phase  I 
program,  was  exercised  to  get  a  preliminary  estimate  of  the  gain  vs.  beam- 
width  tradeoff.  The  code  simulated  the  SNOTEL  transmitters  In  Boise,  Idaho 
and  Ogden,  Utah  and  a  receiver  antenna  array  In  San  Diego  similar  to  that 
proposed  for  the  Phase  II  program. 

The  receiver  antenna  array  will  be  capable  of  producing  three 
beamwldths  of  10°,  16°,  and  40°  with  corresponding  gains  of  18  dB,  16  dB, 
and  12  dB,  respectively.  Note  the  near  exact  Inverse  relationship  between 
beamwldth  and  received  power  (not  gain).  Beamwldths  smaller  than  10°  are 
very  difficult  to  achieve  owing  to  practical  limitations  on  antenna  array 
sizes  and  the  physical  effects  of  the  partially  conducting  ground  plane  of 
the  earth. 

Because  the  code  simulates  random  arrival  of  meteors  and  contains 
numerous  physical  effects,  each  simulation  runs  5  to  8  hours  on  an  IBM  PC. 
The  present  Implementation  of  the  code  on  a  microcomputer  avoids  CPU  costs 
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and  Is  convenient  for  overnight  runs.  Available  time  In  the  Phase  I 
program  did  limit*  however*  the  scope  of  our  predictions.  We  only  had  time 
to  consider  the  receiving  beam  to  be  aimed  directly  at  each  of  the  trans¬ 
mitters  In  turn.  Had  we  aimed  the  receiving  beams  at  the  hot  spots 
Instead*  the  results  would  be  expected  to  look  completely  different. 

The  results  of  the  MBCS  code  are  summarized  In  Table  3.1.  Table  3.1 
shows  the  percent  probability  that  a  meteor  channel  Is  active  In  a  receiver 
at  San  Diego  from  transmitters  In  Boise  and  Ogden  at  0600  and  1200  for 
three  different  ga In/beamwidth  combinations.  The  pairs  of  numbers  In 
parentheses  represent  the  95%  confidence  Interval  about  the  predicted 
probability.  The  last  column  Is  the  probability  of  connectivity  for  two 
antennas  operating  simultaneously  and  looking  10°  to  either  side  of  the 
great-circle  path  between  the  transmitter  and  receiver. 

The  MBCS  results  clearly  indicate  that  when  the  receiver  antenna  Is 
pointed  directly  at  the  transmitter,  wide  low-gain  beams  have  greater 
probability  of  connectivity.  This  result  Is  not  too  surprising,  because 
the  wide  beams  encompass  more  of  the  hot  spots  on  either  side  of  the 
relatively  quiescent  great  circle  path  between  receiver  and  transmitter.  A 
more  Interesting  calculation  that  will  be  done  In  Phase  II  Is  to  focus 
narrow  beams  Individually  on  the  hot  spots.  Then  we  might  find  an  optimal 
beamwldth  perhaps  comparable  to  the  width  of  the  hot  spots. 

The  three  different  galn/beamwldth  combinations  can  be  readily 
achieved  by  selectively  activating  only  certain  antennas  In  each  8-antenna 
array.  This  procedure  was  shown  schematically  In  Figure  2.19. 

The  Yagl  antenna  arrays  are  cheap*  but  unwieldy  and  difficult  to  move 
once  Instal  led.  One  should  plan  to  have  three  8-antenna  arrays  at  the 
receiving  site.  Two  of  the  arrays  would  be  permanently  fixed.  The  third 
array  would  be  rotated  according  to  the  test  plan  (as  Infrequently  as 
possible)  to  perform  the  gain  vs.  beamwldth  experiments  and  angular  and 
spatial  diversity  experiments. 

The  data  from  the  gain  vs.  beamwldth  experiments  should  be  compared 
with  the  MBCS  code  predictions  both  to  validate  the  code  and  to  develop  a 
good  understanding  of  the  gain  vs.  beamwldth  tradeoffs.  The  benefits  will 
be  the  capability  to  design  any  meteor-burst  commun Ications  system  with 
optimal  beamwldth. 
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TABLE  3.1.  PERCENT  PROBABILITY  OF  CHANNEL  CONNECTIVITY 
(AND  95*  CONFIENCE  INTERVAL) 
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3.2  Angular  Diversity 

If  neighboring  antenna  arrays  directed  towards  different  parts  of  the 
sky  could  receive  signals  uncorrelated  with  each  other#  then  this  would 
have  significant  tactical  Implications  for  deployment  of  a  meteor  burst 
communications  system.  It  would  suggest#  for  example#  that  receiving 
arrays  could  be  colocated  and  Increase  the  combined  percent  copy  by 
receiving  Independent  paths. 

The  COMET  system  operating  In  Europe  has  already  shown  that  wlde-beam 
receiving  antennas#  when  directed  at  a  common  volume#  do  not  obtain 
Independent  paths#  even  when  separated  by  35  km.  This  result  does  not 
address  the  question#  however#  of  whether  narrow-beam  antennas,  focused  on 
separate  hot  spots  can  obtain  Independent  paths.  We  exercised  the  Meteor 
Burst  Communications  Simulation  code  to  get  a  preliminary  answer. 

The  results  of  the  MBCS  code  were  surprising#  but  made  sense  In 
retrospect.  Since  the  two  hot  spots  are  typically  separated  by  the  order 
of  100  km,  and  since  typical  meteor  trails  are  tens  of  km  long,  we 
originally  expected  that  signals  received  by  each  of  two  beams  focused  on 
separate  hot  spots  would  be  uncorrelated.  The  MBCS  code  suggests 
otherwise. 

The  code  was  run  again  to  simulate  the  SNOTEL  transmitters  in  Boise 
and  Ogden  and  an  8-antenna  Yagl  array  In  San  Diego.  In  the  simulation, 
four  antennas  In  the  array  produced  a  16°  receiving  beam  at  16  dB  directed 
10°  to  the  left  of  the  great  circle  path  to  the  transmitter#  while  four 
other  antennas  produced  a  similar  beam  10°  to  the  right.  The  main  lobes  of 
the  two  beams  did  not  overlap  within  the  3  dB  points.  The  simulation  was 
done  first  for  the  transmitter  at  Boise  and  then  for  Ogden.  Boise  and 
Ogden  subtend  an  angle  of  20°  from  San  Diego. 

Table  3.2  shows  the  results  of  this  simulation.  P(L)  and  P(R)  are  the 
probabilities  that  the  left  channel  (10°  to  the  left  of  the  great  circle 
path)  and  right  channel  are  active#  respectively.  P(L/R)  Is  the 
probability  that  the  left  channel  Is  active  at  any  time  that  the  right 
channel  Is  known  to  be  active.  P(R/l)  Is  the  reverse.  The  correlation 
coefficient  Is  a  measure  of  the  correlation  of  signals  between  left  and 
right  channels.  In  all  cases#  the  possibility  that  the  left  and  right 
channels  are  uncorrelated  In  less  than  one  part  In  one  thousand. 
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The  reason  that  the  MBCS  code  shows  such  strong  correlations  between 
beams  directed  at  separate  hot  spots  Is  that  the  side  lobes  of  the  beams* 
although  down  by  more  than  about  10  dB  at  the  location  of  the  other  hot 
spot  *  still  have  enough  gain  to  pick  up  many  of  the  meteor  trails  In  the 
other  hot  spot. 

The  MBCS  code,  as  presently  configured  assumes  a  dipole-antenna 
receiving  pattern.  A  more  realistic  receiving  pattern  corresponding  to  the 
Yagl  antenna  array  will  be  put  Into  the  code  In  Phase  II,  but  the  results 
are  expected  to  be  qualitatively  similar. 

The  experimental  configuration  for  testing  these  predictions  was  shown 
In  Figure  2-20.  The  gain  and  beamwldths  of  the  two  beams  will  be  the  same 
as  In  the  simulations.  The  beams  will  be  directed  towards  the  hot  spots  on 
either  side  of  the  transmitter.  Correlations  of  signals  from  different 
transmitters  will  be  observed  as  the  overlapping  of  the  audible  tones  from 
each  of  the  two  receiving  arrays.  Tests  will  be  done  for  transmitters  In 
different  directions  and  at  different  times  of  day  to  Include  the  effects 
of  the  diurnal  cycle  of  the  hot  spot  Intensities. 

Although  the  code  suggests  that  the  signals  from  the  hot  spots  will  be 
correlated.  It  also  suggests  that  there  may  be  an  advantage  of  using  split 
receiver  beams.  The  last  column  of  Table  3.1  showed  that  the  highest 
probability  of  channel  connectivity  Is  achieved  by  using  all  eight  antennas 
In  the  array  to  produce  two  beams  aimed  at  separate  hot  spots.  The 
Increase  In  percent  connectivity  over  the  next  best  configuration,  which 
was  a  wide  beam  directed  along  the  great  circle  path  to  the  transmitter 
ranged  from  5035  to  80%  In  the  cases  examined. 

Perhaps  the  advantage  of  the  split  beam  Is  more  clearly  displayed  In 
terms  of  meteor  bursts  observed  per  hour.  Table  3.3  shows  these  code 
results  for  the  same  simulation  as  In  Table  3.1.  Again,  95%  confidence 
Intervals  are  given.  The  results  show  that  for  a  receiving  beam  directed 
along  the  great  circle  path  to  the  transmitter,  wide,  low-gain  beams  are 
preferable,  but  a  split  beam,  with  each  looking  at  separate  hot  spots, 
picks  up  the  most  meteors  per  hour. 

This  Is  not  to  say  that  a  split  beam  Is  the  optimal  configuration.  By 
no  means  have  all  likely  configurations  been  tried  In  the  simulations.  For 
example,  a  single,  high-gain  narrow  beam  directed  at  a  single  hot  spot  has 
not  yet  been  simulated.  It  Is  quite  conceivable  that  the  optimal 
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TABLE  3.3.  METEOR  BURSTS  OBSERVED  PER  HOUR  (AND  95*  COMMENCE  INTERVAL) 
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configuration  Is  a  single  beam  that  changes  direction  every  twelve  hours  to 
focus  on  the  hotter  hot  spot. 

The  purpose  of  the  angular  diversity  experiments  would  be  to  develop 
an  understanding  of  the  physics  Issues  determining  optimal  antenna 
configurations  at  a  single  receiving  site  through  mutually  validated 
experiments  and  code  simulations. 

3.3  Spatial  Diversity 

The  Issue  of  spatial  diversity  considers  advantages  In  a  receiving 
network  In  terms  of  percent  copy  of  spatially  separated  receiving  antennas. 
Alternatively#  by  reciprocity#  one  may  consider  spatially  separated 
transmitting  antennas. 

Since  the  MBCS  code  suggests  that  colocated  receiving  antennas 
directed  towards  separate  hot  spots  may  pick  up  signals  that  are  stongly 
correlated#  the  next  question  becomes  how  far  apart  must  the  receiving 
antennas  be  In  order  to  pick  up  uncorrelated  signals.  Some  work  on  this 
question  was  done  In  the  COMET  network  In  Europe,  but  In  those  experiments 
the  receivers  all  viewed  a  common  volume.  We  are  more  Interest  In  the 
various  configurations  that  promise  higher  percent  copy  by  selectively 
focusing  on  hot  spots.  For  example#  In  a  two-transmitter  system#  the 
optimal  percent  copy  may  be  achieved  by  directing  one  receiving  antenna 
array  towards  the  left  hot  spot  of  the  left  transmitter#  and  another  to¬ 
wards  the  right  hot  spot  of  the  right  transmitter#  and  depending  on  the 
side  lobes  to  pick  up  the  meteor  bursts  In  the  hot  spots  between.  The 
optimal  receiving  antenna  directions  can  be  expected  to  be  dependent  on 
both  angular  and  spatial  separation  of  the  transmitters#  beamwldth  of  the 
receiving  antennas#  and  time  of  day. 

The  experimental  configuration  to  test  the  correlation  of  signals 
obtained  by  spatially  diverse  receivers  was  shown  In  Figure  2.21.  Three 
receiving  antenna  arrays  would  be  located  at  the  receiving  site.  The  two 
arrays  pointed  at  the  outermost  angular  separations  would  remain  fixed.  A 
third  would  be  rotated  to  point  In  various  directions  between  the  outermost 
two.  This  arrangement  represents  the  fewest  possible  antennas  that  can 
address  the  Issue  of  spatial  diversity  without  requiring  Impractical 
amounts  of  time  for  reorienting  antennas.  Similarly#  the  two  transmitters 
at  the  outermost  angles  would  remain  fixed.  A  third  and  possibly  a  fourth 
transmitter#  which  would  be  operating  simultaneously  with  these  two#  should 


be  located  at  a  variety  of  distances  from  the  receiving  site  and  a  variety 
of  directions  between  the  outermost  sites. 

The  diversity  achieved  In  this  manner  would  enable  a  determination  to 
be  made  of  correlation  footprints  as  functions  of  range  and  angle  as  well 
as  time  of  day  and  beamwldth  and  direction  with  respect  to  hot  spots.  This 
data  Is  critical  to  developing  the  capability  to  optimize  designs  of  meteor 
burst  communications  systems. 

Our  preliminary  designs  of  the  test  systems*  discussed  In  Subsection 
2.3*  were  guided  by  fact  that  the  equipment  for  the  Phase  II  experiments  Is 
simple  and  cheap*  and  that  the  dominating  costs  are  In  labor.  Where 
tradeoffs  were  possible*  we  have  chosen  additional  equipment  to  reduce 
labor  costs  and  overall  costs.  This  Is  reflected*  for  example.  In  the 
redundancy  of  receiving  antenna  arrays. 


4.0  PHASE  II  PROGRAM  PLAN 


4.1  Phase  II  Tasks 

This  subsection  discusses  the  four  principal  tasks  of  the  VHF 
Reception  Program  -  Phase  II. 

TASK  1  -  EXPERIMENT  PREPARATION 

TASK  1.1.  Frequency  A1  location 

The  experiment  requires  the  transmission  and  reception  of  an  RF  signal 
In  the  low  VHF  (30-50  MHz)  range.  It  Is  expected  that  at  any  one  t1me»  no 
more  than  three  transmitters  will  be  operating  with  one  central  receiver 
site.  Two  transmitters  wll  1  be  operating  at  fixed  sites;  the  third  will 
operate  at  one  of  four  or  five  locations.  Frequency  authorizations  will  be 
sought  for  as  many  as  nine  transmitting  sites,  however,  to  allow  for 
difficulties  In  obtaining  some  authorizations  or  site  use  approvals.  All 
sites  will  be  authorized  for  the  same  frequency  band.  A  license  and 
frequency  authorization  to  be  obtained  from  the  FCC  requires  the  submission 
of  documentation  and  technical  Information  Including  specific  locations  and 
characteristics  of  the  transmitters  and  receivers. 

TASK  1.2.  Facilities  Procurement 

The  use  of  government  space  and  facl  1  Itles  requires  the  request  for 
use  and  coordination  of  users  and  hosts  to  satisfy  government  or  military 
requirements  for  support  of  a  government  sponsored  experiment.  Our 
preliminary  discussions  Indicate  that  military  sites  are  willing  to  host 
experiments  provided  that  the  experimenters  supply  the  housing  for 
equipment,  such  as  sheds  or  trailers. 

TASK  1.3.  Test  Planning 

A  detailed  test  plan  designed  to  satisfy  the  test  goals  and  objectives 
of  the  experiment  must  be  developed  which  provides  the  general  direction 
and  specific  Instruction  for  the  fielding,  operation  and  data  evaluation  of 
the  experiment. 

TASK  1.4.  Test  Schedule 

A  test  schedule  will  Identify  the  sequence  of  activities  based  on  the 
definition  of  milestones  and  estimates  of  effort  need  to  accomplish  the 
experiment  tasks. 
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TASK  2  -  EXPERIMENT  FIELDING 

TASK  2.1.  Hardware  Design 

The  hardware  design  will  Include  the  selection  of  commercially  avail¬ 
able  equipment  and  the  Integration  and.  If  necessary,  design  modification 
needed  to  perform  the  experiments  defined  In  the  test  plan. 

TASK  2.2.  Procurement 

Equipment  must  be  procured  based  on  the  equipment  specified  In  the 
hardware  design.  A  preliminary  design  should  be  developed  to  Identify  any 
long-lead  Items.  A  cost  trade-off  should  be  evaluated  to  Identify  the  most 
advantageous  buy-rent  options  based  on  the  hardware  design  and  the  test 
schedule. 

TASK  2.3.  Fabrication  and  Assembly 

While  much  of  the  experiment  hardware  will  be  off-the-shelf 
electronics,  the  Integration  and  hookup  of  the  Interfacing  units  and  the 
design  of  selected  Items  will  require  the  fabrication  and  assembly  of  the 
units  Into  an  operable  system. 

TASK  2.4.  Installation  and  Checkout 

The  Installation  of  the  equipment  for  the  VHF  reception  experiment 
will  require  not  only  the  assembly  and  hookup  of  the  transmitters  and 
receivers  and  antennas,  but  will  include  accurate  surveys  to  establish  the 
accurate  pointing  of  antennas  and  antennas  arrays. 

TASK  3  -  TEST  OPERATIONS 

The  proposed  VHF  reception  experiment  Includes  three  different 
experimental  efforts  that  use  the  antenna  arrays  to  evaluate  different 
experimental  configurations.  The  present  concept  Includes  the  capability 
to  simultaneously  receive  3  different  signal  paths,  and  Includes  the 
options  to  vary  the  galn/beamwldth  of  the  antenna  arrays.  Each  experiment 
will  require  the  specific  configuration  of  transmitters  and  receiver 
antenna  arrays  defined  In  the  test  plan.  Data  will  be  acquired  and 
recorded  on  a  scheduled  basis.  Upon  completion  of  each  experiment  phase, 
the  antennas  will  have  to  be  reconfigured  to  the  orientation  and  configura¬ 
tion  needed  for  the  next  experiment. 

Transmitters  will  be  remotely  operated  (turned  on  and  off)  using 
normal  dial-up  telephone  control  from  the  central  receiver  site. 
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The  specific  test  and  an  evaluation  of  test  results  will  determine  the 
selection  of  options  to  move  one  or  more  of  the  remote  transmitters  to 
different  authorized  and  licensed  locations. 

TASK  4  -  TEST  RESULTS 

An  objective  of  this  proposed  VHF  reception  experiment  Is  to  validate 
existing  meteor  channel  computer  models  and  propagation  theory.  A  set  of 
preliminary  predictions  which  model  the  experiment  configurations  will  be 
compared  to  the  experimental  data  acquired  during  the  test  operation  phase 
of  the  program.  A  continuing  analysis  effort  will  be  maintained  while  test 
operations  are  In  progress.  Quick-look  analysis  will  provide  an  evaluation 
of  the  usefulness  of  the  experimental  data  and  will  be  used  to  make 
adjustments  to  the  test  plan  to  optimize  the  test  configuration  and  data 
format.  An  assessment  of  the  test  results  and  a  comparison  with  predic¬ 
tions  will  be  described. 

4.2  Statement  of  Work 

A  statement  of  work  for  the  Phase  II  program  has  been  developed  as 
follows. 

TASK  1  -  EXPERIMENT  PREPARATION 

Perform  preliminary  authorization  activities  necessary  for  the  VHF 
Reception  Program  -  Phase  II#  Including  preparation  of  documentation  and 
coordination  needed  to  obtain  proper  frequency  authorization  and  government 
licensing  of  radio  transmitters  and  radio  sites#  and  approval  for  the 
Installation  of  transmitters  and  receivers  on  government  or  military 
property#  and  use  of  serv  Ices  and  facl  Titles  needed  for  test  operations. 
Develop  a  detailed  test  plan  and  schedule  for  the  performance  of  the  VHF 
reception  experiment.  The  plan  will  provide  for  the  sequencing  of  the 
experiment  parts  to  ensure  that  cyclical  effects#  such  as  diurnal  and 
seasonal  variations#  are  Included  in  the  acquired  data  base. 

TASK  2  -  EXPERIMENT  FIELDING 

Perform  the  fielding  of  the  experiment#  Including  the  design#  procure¬ 
ment#  fabrication  and  assembly  of  the  required  hardware  and  electronics. 
The  equipment  will  be  Installed  and  checked  out  In  accordance  with  the  test 
plan  and  schedule. 


TASK  3  -  TEST  OPERATIONS 

The  test  operations  will  be  performed  In  accordance  with  the  test 
plan.  Operations  wl 1 1  Incl ude  control  of  the  remote  transmitter  sites, 
receiver  operation,  and  data  recording.  The  test  facilities  will  be  recon¬ 
figured  as  scheduled  for  each  of  the  various  experiments. 

TASK  4  -  TEST  RESULTS 

Test  results  will  be  evaluated  by  comparing  preliminary  modeling  and 
prediction  efforts  to  test  data.  Assessment  of  these  data  will  be 
documented  In  reports. 

4.3  Schedule 

A  detailed  test  plan  which  defines  exact  test  configurations:  path 
orientation  (surveyed  angular  measure  from  true  North),  antenna 
combinations,  transmit  on-off  times,  and  data  collection  criteria  will  be 
developed  based  on  final  site  locations  and  analytical  requirements.  A 
representative  schedule  showing  estimated  hardware  design,  procurement  and 
fabrication  times,  allocation  of  antenna  reconfiguration  times  and  a 
sequence  of  dally  and  seasonal  experiment  data-taklng  periods  Is  shown  In 
Figure  4.1. 

4.4  Manpower  Allocation 

Estimates  of  man-months  that  should  be  allocated  to  each  activity  In 
Phase  II  are  displayed  In  Table  4.1. 

TABLE  4.1.  MAN-MONTHS  ALLOCATED  TO  PROGRAM  ACTIVITIES 


ACTIVITY  MAN-MONTHS 

Program  Management  4 
Analysis  6 
Experiment  Design  4 
Installation  and  Maintenance  8 
Test  Operations  20 
Data  Reduction  8 
Site  Selection,  Coordination  3 


TOTAL 


53 


Ffgure  4.1.  Schedule  for  VHF  Reception  Program  —  Phase  II 
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